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Abstract The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful 
tool for observing the Sun at high spatial, temporal, and spectral resolution. These capa¬ 
bilities can address a broad range of fundamental scientific questions in solar physics. The 
radiation observed by ALMA originates mostly from the chromosphere - a complex and dy¬ 
namic region between the photosphere and corona, which plays a crucial role in the transport 
of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. 
Based on first solar test observations, strategies for regular solar campaigns are currently be¬ 
ing developed. State-of-the-art numerical simulations of the solar atmosphere and modeling 
of instrumental effects can help constrain and optimize future observing modes for ALMA. 
Here we present a short technical description of ALMA and an overview of past efforts and 
future possibilities for solar observations at submillimeter and millimeter wavelengths. In 
addition, selected numerical simulations and observations at other wavelengths demonstrate 
ALMA’s scientific potential for studying the Sun for a large range of science cases. 

Keywords Sun • photosphere • chromosphere • corona • magnetohydrodynamics • radiative 
transfer • flares • prominences 


1 Introduction 

Despite decades of intensive research, the chromosphere of our Sun - an important atmo¬ 
spheric layer between the photosphere and the corona - is still elusive owing to the com¬ 
plicated nature of the partially ionized chromospheric plasma and its intricate interaction 
with radiation and the arising technical challenges for observing the emergent radiation (see 


C. E. Alissandrakis 

Department of Physics, University of loannina, GR-45110 loannina, Greece. 

A. Nindos 

Department of Physics, University of loannina, GR-45110 loannina, Greece. 

S. K. Solanki 

Max-Planck-Institut for Sonnensystemforschung, Justus-von-Liebig-Weg 3, 37077 Gottingen, Germany 
School of Space Research, Kyung Hee University, Yongin, Gyeonggi, Republic of Korea 

L. Rouppe van der Voort 

Institute of Theoretical Astrophysics, University of Oslo, Postboks 1029 Blindem, N-0315 Oslo, Norway 

M. Shimojo 

National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 
Y. Kato 

Institute of Theoretical Astrophysics, University of Oslo, Postboks 1029 Blindem, N-0315 Oslo, Norway 

T. Zaqarashvili 

Institute of Physics, University of Graz, UniversitLtsplatz 5, 8010 Graz, Austria 

Abastumani Astrophysical Observatory at Ilia State University, University St. 2, Tbilisi, Georgia 

E. Perez 

SUPA, School of Physics & Astronomy, University of Glasgow, UK 
C. L. Selhorst 

University of Vale do Paraiba (UNIVAP), Av. Shishima Hifumi, 2911 - Urbanova - CEP 12244-000, Sao Jose 
dos Campos, Sao Paulo, Brazil 

M. Barta 

European ARC, Czech node. Astronomical Institute ASCR, Ondrejov, Czech Republic 
Astronomical Institute, Academy of Sciences, Ericova 298, 251 65 Ondrejov, Czech Republic 



4 


06^2014: AR12157 (x^)-(-73r.-313') 274ui!-2010: AR11084, (x^)*(-S2ir,-33(r) 



24-Jun-2014: coronal bole, (x,y)>('129*,-8r) 14.Jufi-2nOH : ARI099fl. (x,y)3(.13ir,-159*) 



X [AfOMC) X [«raee] 


Fig. 1 Examples of chromospheric images taken in the core of the Ca II infrared triplet line at A = 854.2 nm 
with the CRISP instrument at the Swedish 1-m Solar Telescope (SST). (a) An active region with an on-going 
Ml.l class flare, (b) The chromosphere above a sunspot close to the limb, (c) A quiet region inside a coronal 
hole close to disk center, (d) A decaying active region with a pronounced magnetic network cell close to disk 
center. 


Fig.[T]l. While advances in instrumentation during recent years have led to improved chromo¬ 
spheric diagnostics, sparking renewed attention for this important interface layer, the diffi¬ 
culties with interpreting the few currently accessible diagnostics with complicated formation 
mechanisms have hampered progress. In contrast, radiation at millimeter wavelengths has 
been known to have great diagnostic potential but could not be exploited to full extent due to 
technical limitations so far, in particular, because technically achievable telescope apertures 
had been small for the relatively long wavelengths. Consequently, single dish observations 
did not have sufficient temporal or spatial resolution, and small interferometric arrays could 
not measure the full range of spatial scales in the chromosphere. 

The Atacama Large Millimeter/submillimeter Array (ALMA), which recently com¬ 
menced operation, is now an enormous leap forward in terms of spatial resolution at mil- 
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Fig. 2 Reconstructed interferometric images of a Quiet Sun region (left) and active region (right) from ob¬ 
servations with the Berkeley-Illinois-Maryland Array (BIMA). The yellow circles in the lower left comer of 
each panel represent the restored beam with a diameter 10”. Adopted from|Loukitcheva et al.|(2009]|2014). 


limeter and submillimeter wavelengths, offering an unprecedented view of the solar chromo¬ 
sphere and its interface with the upper photosphere and low transition region. For extended 
sources like the Sun, ALMA will achieve a spatial resolution that is close to what is currently 
possible at visible and infrared wavelengths. In addition, ALMA observes with high tempo¬ 
ral resolution and can change between different receivers and thus wavelength bands, which 
are subdivided in a large number of spectral channels, both polarised and unpolarised. The 
major advantage of ALMA compared to other diagnostics is that it serves as a nearly linear 
thermometer of the solar plasma because the observed brightness temperature is closely re¬ 
lated to the gas temperature of the continuum forming layer. It thus probes the local plasma 
conditions that are, at least in principle, much easier to interpret. Hence, ALMA will be 
strongly complementary to other ground-based and space-borne observatories. Coordinated 
observations with other observatories will therefore provide new science opportunities to 
improve our understanding of the dynamics and energetics of this highly complex interface 
between the solar photosphere and the corona. 

ALMA’s unprecedented capabilities will allow addressing a wide range of topics and 
major questions in contemporary solar physics ( |Bastian| |2002[ [Karlicky et al.|[20TTj |Benz| 
|et al.|2012a[|Wedemeyer et al.|2015| l, including the dynamics, thermal structure and energy 
transport in the “quiet” solar chromosphere, active regions and sunspots, spicules, promi¬ 
nences and filaments, and flares. An example of the challenging complex structure of the 
solar chromosphere is shown in Fig. [T] These high-resolution images were obtained by the 
CRISP (CRisp Imaging SpectroPolarimeter) instrument at the Swedish I-m Solar Telescope 
(SSlQl in the core of the Call infrared triplet line at a wavelength of 854.2 nm. This spec¬ 
tral line is among the currently most used chromospheric diagnostics. The images exhibit 
structures on spatial scales down to the resolution limit of the telescope (~0.”2 at 854.2 nm 
but down to ~0.”1 for, e.g., Ca II K). In the millimeter range, however, the spatial resolu¬ 
tion of observations has been significantly lower in the pre-ALMA era as can be seen from 
the otherwise elaborate examples in Fig. which have been obtained with the Berkeley- 
Illinois-Maryland Array (BIMA) at a wavelength of T = 3.5 mm ( [Loukitcheva et al.|2009[ 

* Please note that the abbreviation SST is used for both the Swedish 1-m Solar Telescope and the Solar 
Submillimeter Telescope. 
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Fig. 3 Full-disk maps of the Sun taken on September 7th, 2014. Left: Brightness temperature map obtained 
by scanning the solar disk with a single ALMA total power antenna (PMOl) at a frequency of 230 GHz 
(/i = 1.3 mm) in a double-circle pattern, which was completed in roughly 10 rain. Right: Corresponding Her 
map. Courtesy of |Phillips et al.|(201 5^. Please note that this commissioning data is shown here for illustration 
purposes only in accordance with ALMA data policy. 


|2014[ l. The magnetic network is clearly visible in the image of a quiet Sun region (left panel) 
but no features smaller than the restored beam size of 10” can be discerned. Consequently, 
the intricate small-scale chromospheric structure inside the network regions, which can be 
clearly seen in the Call image in Fig.[^, remains unresolved in the BIMA maps. This is 
true for all other instruments operating at millimeter wavelengths, so that much of the rich 
chromospheric fine-structure could not be observed at these long wavelengths so far. The 
situation now changes completely with ALMA, which opens up a new domain of spatial 
resolution at millimeter wavelengths. 

In preparation of regular solar observations, which are expected to start in Cycle 4 in 
2016 after the initial solar observing modes have been demonstrated (see Fig.[^, this review 
article provides a summary of ALMA’s capabilities and outlines potential science cases. An 
introduction to the origins of the radiation at millimeter and submillimeter wavelengths in 
the solar atmosphere is given in Sect. followed by a description of the technical parame¬ 
ters and capabilities of ALMA and its diagnostic potential for solar observations in Sect.[^ 
Sectionj^provides a broader and more detailed overview of potential science cases. Conclu¬ 
sions and an outlook are given in Sect.[^ This article was produced as a joint activity of the 
Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMOlsj^ 
http://ssalmon.uio.no), which is an international network aiming at anticipating key solar 
science with ALMA through simulation studies l|Wedemeyer et al.|2015|l. 


^ The network is open to everybody who has a professional interest in contributing to potential ALMA so¬ 
lar science, which include or require simulations. More information and registration at http://ssalmon.uio.no 
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2 The formation of millimeter radiation in the Sun 


In the absence of strong magnetic fields, e.g. in quiet Sun regions, the opacity at millime¬ 
ter wavelengths is mainly due to electron-ion free-free absorptioip] and - to lesser extent 


- H free-free absorption (cf., e.g.. White et al. 2006 Shibasaki et al. 2011 and refer¬ 


ences therein). The resulting thermal radiation and their diagnostic potential are described 
in Sect. |2.1| Radio recombination lines and molecular lines at (sub-)mm wavelengths are 
discussed in Sects. |2.2| and |2.3| respectively. As detailed in Sect. |2.4| the gyroresonance 
process does not contribute significantly to the thermal emission at millimeter wavelengths, 
while the gyrosynchrotron mechanism can produce non-thermal emission in connection with 
magnetic fields and high-energy particles. Additional emission mechanisms, which are dis¬ 
cussed in connection with a yet unexplained high-frequency (i.e., sub-THz) flare emission 
component (see Sect. |4.4| and |Fleishman & Kontar|2010[|Krucker et al.|2013^ include syn¬ 
chrotron emission from relativistic positrons and electrons, diffusive radiation, and possibly 
Cherenkov emission hut require further investigation. 


2.1 Thermal continuum radiation at millimeter and submillimeter wavelengths 


The two major opacity sources relevant for thermal continuum radiation at millimeter and 
suhmillimeter wavelengths in the solar atmosphere are due to free-free absorption of the 
following two kinds: (1) Electron-ion free-free absorption^ (H"^ + e~ + y —> -r e“*, i.e. 

interaction of a free electron with radiation in the Coulomb field of a charged ion). (2) H“ 
free-free absorption (ff -r e“ -H y —> H -r e“*), occurring in the field of a neutral hydro¬ 
gen atom which is polarized by the Coulomb field of the passing electron. At millimeter 
wavelengths, free-free absorption can be treated either by quantum as well as - by virtue 
of the correspondence principle - classical theory (see |Oster|1961a[ l. Due to its importance 
a lot of work was dedicated to the treatment of this process. Different approaches and ap¬ 
proximations led to a multitude of prescriptions of the absorption coefficient whose ranges 
of validity and accuracy are not always obvious to the end-user (for a discussion, see |Oster] 
|1961b| l. However, the possibility to treat the problem fairly rigorously leaves these ambigu¬ 
ities on a level which is typically immaterial for application in solar physics. We point the 
reader to the widely used quantum mechanical calculations of the Gaunt factor by |Karzas &| 
|Latter| ( [l961| l. 

Here, we follow [DuIk| ( |1985l l and write the free-free absorption coefficient due to electron- 
ion free-free absorption^ in the absence of magnetic fieltj^semi-classically as 

1 (lyyl 4;re'* V ■ 72 „ ^ 

Aions.ff ~ a..! — I '-r\^n / j ^^.?ions,tf(T, v) (1) 


v2 


where gions.ff is the thermally averaged Gaunt factor]^ 


^ V3, / 

gions,ff(T, y) = -In ———^ 


(2) 


^ This process is also referi'ed to as “inverse thermal bremsstrahlung” in the sense that it is the absorption 
process reverse to photon emission by bremsstrahlung (e.g., |Frolov|2Q15| . 


^ See Sect. 


2.5 


for the magnetized case. 


Note, that Eq. 18 of 


Dulk 

1985 

18 as 

is. It is 


is likely in error, the numerical value of 18.2 of Bulk’s Eq. 20 cannot 


be reproduced using Eq. 18 as is. It is likely that the Gaunt factor should read as given in Eq.[^here which 
brings Bulk’s Gaunt factor much closer to results of |Karzas & LatterH1961) and |Oster|jl970| . The numerical 
evaluation given later is corrected accordingly. 
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Fig. 4 Opacity at 1 mm wavelength due to electron-ion free-free absorption and free-free absorption as 
a function of temperature in units of the Thomson cross section for electron scattering = 6.65 X 

10“^^ m^). The opacity was normalised to the electron density and to the effective density of absorbing 
ions or atoms, i.e. to Z^n, in the case of electron-ion bremsstrahlung and to the number density of neutral 
H-atoms in the case of . Please note that the opacity is affected by the low hydrogen ionisation and 
thus low election density at low temperatures. 


valid for temperatures T < 2 x 10^ K. c is the speed of light, F x 1.781 is Euler’s constant, 
is the Boltzmann constant, v the photon frequency, T the (kinetic) electron temperature. 
We the electron density, n, the density of ion i of charge Z,, e the electron charge (4.8... x 
10“'° esu), nie the electron mass, and 



(3) 


the plasma frequency. Numerically, Eq. [^evaluates to (in cgs units) 

Elions,ff ^ 9.78 X 10-3-^ X (17.9 + InZ^/^ _ iny) 


[cm 


which leads to the following simplified expression ( |Kundu|1965[ l 

r? 

X = ^iT,ey 


1 / 273/2 


(4) 


(5) 


where = tii is assumed, n is the index of refraction and e) is a slowly varying function 
of electron temperature and density, which evaluates to - 0.12 in the chromosphere and 
- 0.2 in the corona. 

Eigure]^ shows the resulting opacity for a reference wavelength of 1 mm together with 
the free-free opacity due to H“ according to |Stallco^ l [1974| l. In a plasma mainly composed 
of hydrogen, electron-ion free-free absorption-* is the main opacity source in the considered 
temperature range as long as the ionisation fraction of hydrogen is not very small. Chromo¬ 
spheric electron densities are in the range «e ~ 10° ... 10" cm“^, implying densities of ions 
and neutral hydrogen of similar order. As a consequence, Thomson scattering - as opposed 
to the solar corona - is an unimportant opacity source in the chromosphere. 
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According to Eq.[^the plasma frequency increases with electron density. In a region of 
rather high electron density of 10" cm“^ the plasma frequency corresponds to a wavelength 
of ~ 10 cm. For a wave of a rather long (by ALMA standards) wavelength of 1 cm this 
results in an 0.5 % reduction of the index of refraction from one. This implies that to good 
precision light rays at ALMA wavelengths travel on straight lines. Moreover, the effect on 
the absorption coefficients is minor since it scales inversely proportional to the index of 
refraction. 

Since the opacity is due to particle interactions, the assumption of LTE conditions holds 
and source function of the radiation continuum is Planckian . At millimeter and radio wave¬ 
lengths A (and for chromospheric temperatures T), the Rayleigh-Jeans limit of the Planck 
function can be used in good approximation (since he!A <K ke T) so that the spectral ra¬ 
diance, i.e. the energy emitted per unit area of emitting surface, per unit time, per unit fre¬ 
quency and per unit solid angle, due to thermal emission is given by 


5,1 (T) = 2hc^ 


e T' 



IcksT A-* 


( 6 ) 


or 


ByiT) 



( 7 ) 


The emergent specific intensity ly is then derived from integration of the source function 
(here the Planck function {T{t)) or By over the optical depth r along the line of 

sight: 

(v) 

Iy= By(T(f))e-Uf, (8) 

Jo 


where (i') is the maximum of the optical depth range at the considered frequency (or 
wavelength) where significant contributions to the source function occur. In the optically 
thick case Tmax = 

In the radio and millimeter range, it is more common to use the brightness temperature 
Tb instead of the continuum intensity ly or I,i, which are related in the following way (as 
derived from Eqs. m- 

Tb = :rr- Iv = :A-r- h. ( 9 ) 


c 


2 c kn 


A common unit for the flux density at these wavelengths is Jansky with 
1 Jy = 10“^® Wm“^ Hz“* = 10“^^ erg s“*cm“^Hz“*. In solar physics, however, it is more 
common to use solar flux units (1 sfu= 10''Jy). Typical values for the Quiet Sun flux den¬ 
sity integrated over the solar disk are on the order of ~ 1.2 10® sfu at 0.3 mm (1 THz), 
~ 1.2 10® sfu at l.Omm (300GHz), ~ 1.5 lO"* sfu at 3.0mm (lOOGHz), and ~ 2 10® sfu 
at 9.0 mm (33 GHz), respectively, based on observed brightness temperatures compiled by 
|Loukitcheva et al.| ( |2004[ l. The flux from active regions (the “S-component”) can be compa¬ 
rable, whereas flares are much brighter. The extreme case shown in Fig. [T^ even reaches a 
peak value of ~ 5 lO"* sfu at 3 mm (100 GHz). 

Numerical modelling of the solar atmosphere and the emergent radiation has proved 
useful for predicting which properties of the radiation can be exploited as plasma diag¬ 
nostic, how to set up the instrument, and how to interpret the resulting data. For instance, 
radiative transfer calculations have been carried out based on model atmospheres for dif¬ 
ferent regions, including the quiet Sun, coronal holes, filaments and active regions ( |Brajsa| 
|et al.|2007||2009[[Romstajn et al.|2009} . The resulting synthetic brightness temperatures for 
a wavelength of 8 mm are in good agreement with corresponding observational results. A 
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Fig. 5 Contribution functions of the continuum intensity at four selected wavelengths between 0.3 mm and 
9.0 mm at solar disk-center based on a numerical model snapshot of the solar atmosphere under quiet Sun 
conditions. For this model (see model E described in Sect. |4.2.'^ , time-dependent hydrogen ionisation and 
with it non-equilibrium electron densities have been taken into account. The contribution functions have 
been averaged horizontally and normalised to their maximum. The dashed vertical line labelled “classical 
temperature minimum” marks the boundary between the photosphere and the model chromosphere. Data 
adapted from|Wedemeyer-Bbhm et ^j2007). 


more extensive study for a broader wavelength range from 0.3 mm to 10.0 mm is currently 
in progress l |Brajsa & et al.|20T5l l. 

The linear dependence on the gas temperature means that ALMA effectively acts as a lin¬ 
ear thermometer of the solar chromosphere and provides measurements of localised plasma 
conditions, which are in principle much easier to interpret than other chromospheric diag¬ 
nostics in the ultraviolet, visible or infrared spectral range. However, the opacity depends on 
the number of free electrons and protons and thus on the ionisation degree of hydrogen. The 
millimeter continuum radiation is therefore affected by deviations from the ionisation equi¬ 
librium of hydrogen due to hnite recombination timescales in the chromosphere. Departures 
from Maxwellian electron velocity distributions may also have detectable effects (see, e.g., 
Fleishman & Kuznetsov|[2014[ l, including not only flare-related kappa distributions (e.g.. 
Oka et al.|2015 I, but also equilibrium distributions (e.g.,|Bian et al.|2014^. 


Furthermore, the opacity increases with wavelength so that the atmosphere becomes 
optically thick at increasing height. In other words, the effective formation height of the 
millimeter continuum radiation increases with height from the temperature minimum at the 
shortest wavelengths to the high chromosphere (and transition region) at the longest wave¬ 
lengths accessible with ALMA ( |Wedemeyer-B6hm et al.|2007| l. This effect is clearly visible 
in the results of the brightness temperature synthesis based on 3D numerical simulations of 
quiet Sun regions, which are described in more detail in Fig. |4.2.2| 

The contribution functions of the continuum intensity at ALMA wavelengths is shown 
in Fig.j^for a selected model, for which important non-equilibrium effects have been taken 
into account (see model E described in Sect. |4.2.2 for more details). These contribution 
functions (published originally by |Wedemeyer-B6hm et al.|2007[ l describe how much of the 
emergent radiation is contributed over height (here for a view directly from the top, thus cor¬ 
responding to solar disk centre). On average, the contribution maximum occurs at a height 
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of 490 km for A = 0.3 mm, 730 km for A = 1.0 mm, 960 km for A = 3.0 mm, and 1170 km 
for A = 9.0 mm in their quiet Sun model, respectively. For comparison, [Loukitcheva et al.] 
( |2015a[ l find average formation heights of 700 km for A = 0.4 mm, 1000 km for 4 = 1 mm, 
1550km for A = 3mm, 1700km for A = 4.5mm and 2000km for A = 9mm, respec¬ 
tively, based on a model of an enhanced network region (see Sect. |4.2.^ . The height of the 
maximum contribution (zm) varies with spatial position and view angle due to the substan¬ 
tial thermal inhomogeneities in the chromosphere. In the quiet Sun model by |Wedemeyer-| 
|Bohm et al.| ( |2007| l, the standard deviation Zm{x,y) ranges from ~ 90 km for A = 0.3 mm 
up to ~ 160 km for A = 9 mm, which is ^Zm(4)/Zm,avg(d) ~ 15 % of the average height 
of the contribution maximum for all wavelengths. Alternatively, the FWHM of the average 
contribution functions are 220 km, 360 km, 440 km, and 480 km, respectively. 

It is important to note that the formation height can thus vary substantially in space 
but also in time ( [Loukitcheva et al.| 2004[|Wedemeyer-Bdh m et al.|2007}|Loukitcheva et ah] 
|2015a^ on time scales of a few 10 s when the variation is due to propagating shock waves. 
On the other hand, accounting for the time-dependent non-equilibrium behaviour of hy¬ 
drogen ionisation leads to less extended formation height ranges compared to the local 
thermodynamic equilibrium (LTE) conditions (see Sect. |4.2.2| for more details on small- 
scale fluctuations). In summary, the continuum intensity from quiet Sun regions at ALMA 
wavelengths originates from relatively narrow layers, whose effective heights increase with 
wavelength, although with a significant overlap in the contribution functions. This property 
allows ALMA to observe different chromospheric layers by changing the wavelength. 

In the general magnetized case, the brightness temperature spectrum noticeably re¬ 
sponds to the magnetic field along the line of sight in the continuum forming layer due 
to the dependence of the free-free opacity on the local magnetic field strength (see Sect. |2.5| 
for details). This phenomenon offers a straightforward way to measure magnetic fields based 
on the observed degree of polarisation in the spatially resolved spectrum. 


2.2 Recombination lines 


A transition between Rydberg states of an atomic ion yields a line at radio frequencies 


V 


2 c Rm 


An 


( 10 ) 


for an ion of effective charge Z (Z = 1 for H I, 2 for He II, etc.) and neutral atomic mass M 
in atomic mass units (amu), where c is the speed of light, n is the principal quantum number 
of the lower state of transition (n -H An) —> n, and Rm is related to the Rydberg constant 
Roo = 109737.315685 cm-' by 


Rm - 


Roi 


1 -H (1822.88848 M-Z)- 


( 11 ) 


Rydberg states of atoms are populated by recombination processes in plasmas and such 
recombination lines have been widely observed in photoionized nebulae and interstellar 
matter. jPupreej ( |T968[ l proposed that such lines might be seen in the solar corona where di- 
electronic recombination naturally overpopulates excited states of abundant ions like Fe XV 
and O VI. Searches for coronal lines at radio frequencies have remained fruitless ( [Bergerj 
|& Simon|1972| l, the main reason being that Rydberg transitions with « > 100 suffer severe 
collisional line broadening that merges them with the continuum even at coronal densities 
( jPolginov & Yakovlev|1974]|Greve|1977[|Hoang-Binh|1982[|Hoang-Binh et al.|1987| l. 
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Meanwhile, limb-brightened, Zeeman-split, mid-infrared emission lines at wavelengths 
near 12/rm ( |Murcray et al.|1981[ |Brault & Noyes|1983[l were identified as Rydberg transi¬ 
tions « = 7 - 6 of Mg I and A11 arising in the low chromosphere ( [Chang & Noyes|1983^ 
and subsequently modelled in some detail i jCarlsson et al.||T992[ [Bruls et al.||1995| l Bor-j 
jeiko & Clarkj ( [T986l l detected emission lines of HI (n = 13,14,15) in the far-infrared, 

V = v/c ~ 25 - 90 cm“’, with a balloon-borne Fourier-transform spectrometer (FTS). Fi¬ 
nally, sub-mm-wave measurements with a FTS at the James Clerk Maxwell Telescope (JCMT) 
demonstrated that recombination lines of HI are detectable near the solar limb at frequen¬ 
cies accessible to ALMA: n = 20-19 at 888.047 GHz ( [Clark et al.|2000a| l and n = 22-21 at 
662.404 GHz ( [Clark et al.[2000b| l. The 21a line appeared at a level of approximately 11 % of 
the continuum with a width ten times broader than thermal Doppler broadening expected at 
the temperature of the lower chromosphere. These pioneering sub-millimeter observations 
had limited angular resolution (19" at 662 GHz). Observations with ALMA will provide 
much better sensitivity, spectral and spatial resolution, and fidelity of line profiles. 

The brightnesses of spectral lines and continuum will provide direct probes of the tem¬ 
perature structure in the chromosphere and transition zone while the line profiles can be 
used to trace microturbulence, flows, and density. Higher in the chromosphere, lines of He I 
might become detectable. With increasing angular distance from the limb, it might even be 
possible to measure the long-sought coronal lines of highly ionized species. In this connec¬ 
tion it is worth noting that coronal O VI 87a at 353.5874 GHz could be observed in the same 
tuning as chromospheric HI 26a at 353.6228GHz, and many other such combinations are 
possible. 

Observations of Rydberg transitions will provide exquisite tests of models of the upper 
photosphere, chromosphere, and corona. Detailed comparisons with models will also de¬ 
mand improved understanding of interesting physics, especially regarding collisional line¬ 
broadening, which remains controversial for low-density plasmas (e.g., [Peac hl|20T4||Hi71 
[M4l i. 


There are other types of atomic transitions that occur at sub-mm wavelengths. Consider, 
for example, the forbidden fine-structure transition in [C III] ls^2s2p - ^P° at 710.2 GHz 
frequency. The upper state of this transition is also the upper state of the CIII 190.9 nm 
line which appears in the chromospheric ultraviolet spectrum of the Sun and solar-type 
stars. The fractional abundance of maximizes at a temperature T = 10^ K. In the 
solar atmospheric model of [Fontenla et al.[ ( [1993| l, this occurs at a height h ~ 2210 km 
above the visible photosphere, where the electron density is approximately 9 x 10®cm“^. 
The excitation of the five lowest levels of is calculated for these conditions, with use 
of the collision strengths and transition probabilities of [Fernandez-Menchero et al. [([2014b. 
Such a simple model that reproduces the observed intensity of the C III 97.7 nm line ( |Warren[ 
[et al.|200T] l implies that the 710 GHz line should appear weakly in absorption against the 
sub-mm-wave continuum that arises deeper in the atmosphere, with an estimated flux of 
-2.3 mjy when averaged over a 5 arcsec diameter on the solar disk. Finally, the positronium 
analog of the 21-cm line lies at 203 GHz and may appear under flare conditions (e.g., [Ellis| 
|& Bland-Hawthorn[2009[ l (see the end of Sect. [4.4. Fj l. 


2.3 Molecular lines of carbon monoxide 

As further detailed in Sect. [4.2. 1[ carbon monoxide molecules have been found to exist in 
cool pockets in the upper photosphere and low chromosphere, both of which are dynami¬ 
cally changing due to the interaction of propagating shock waves and magnetic fields. CO 
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PH0ENIX3O spectrum frx ALMA band 6 
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Fig. 6 Synthetic spectrum for the wavelength range covered by ALMA band 6 (left) and a close-up for a 
spectral line at a wavelength of d = 1.3908 mm (right). The spectrum has been calculated in full 3D with 
the radiative transfer code PHOENIX/3D jHauschildt & Baron|2010) in LTE but includes scattering. For 
this preliminary result, 33 X 33 columns from the 3D magnetohydrodynamic model by jWedemeyer-Bohm] 
|et al.|20i^ has been used as input. The grey area marks the range in brightness temperature covered by all 
considered positions in the model, whereas the thick black line represents the horizontally averaged spectrum. 
The thick dashed lines mark the average plus / minus the standard deviation. The right panel contains a few 
exemplary line profiles, too. 


might be directly detectable with ALMA because rotational transitions of the CO molecule 
lie throughout the millimeter and sub-millimeter parts of the spectrum. At 3500 K, the CO 
J = 6 —> 5 and 1 —> 6 lines at 691 GHz and 806 GHz, respectively, might have brightness 
temperatures of the order of 20 K, to be consistent with the observed equivalent widths and 
line widths of vibration-rotation lines analyzed by |Ayres et al.| l [2006| l. The line widths are 
expected to be large due to pressure broadening. Thus calibration and baseline subtraction 
will be a challenge. However, we know today that the chromosphere exhibits substantial 
temperature fluctuations on short time-scales and on small spatial scales, resulting in corre¬ 
sponding fluctuations of the CO number densities as well as the continuum optical depths. 
It is difficult to estimate to which extent CO rotational lines will be detectable with ALMA 
in the presence of such fluctuations. Numerical simulations involving CO (e.g., Asensio| 


[Ramos et al.|2003j|Wedemeyer-B6hm et al.|2005[|200^|Wedemeyer-B6hm & Steffm 2007^ 

will help answer this question. Such detailed predictions of CO line profiles would be key 
to prepare and interpret ALMA observations. Observations of molecular lines of CO are 
interesting as regarding the thermal structure of the solar chromosphere, but also for the 
determination of elemental abundances and isotopic ratios (see, e.g., [Ayres et al.|2013[ and 
references therein). 


A first attempt at predicting detailed ALMA spectra is presented in Fig.l^ The 3D mag¬ 
netohydrodynamic model by [Wedemeyer-Bohm et al.[ P012[ see Sect. |4.2^[ i nch Table [T 
andFig.|10[l has been used as input for the radiative transfer code PHOENIX/3D l [Hauschildt 
[& Baron[[2010[ [Berkner et al.[[2013^ . The resulting radiative flux for each of the selected 
33 X 33 model columns for 30,000 wavelength points was converted to brightness tempera- 
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tures for Fig.[^ Already this early result shows that the continuum level is expected to vary 
significantly whereas the included spectral lines are very narrow. The line depths are on av¬ 
erage up to ~ 120K and in extreme cases up to 250 K. The calculations have so far been 
carried out under the assumption of LTE for the occupation numbers and the molecular CO 
number densities, whereas background continuum scattering is included. 

The LTE assumption for the occupation numbers is more likely to be valid, owing to 
the fast collision rates (with pervasive atomic hydrogen) that rapidly thermalize the pure 
rotational levels, especially given the long radiative decay times. The instantaneous chemi¬ 
cal equilibrium assumption for the CO number densities is more questionable, owing to the 
increasing difficulty of making the molecules at progressively lower densities at higher alti¬ 
tudes above the photosphere; and the likely frequent interruptions of the formation process 
due to passing hot waves and shocks. Future detailed non-LTE calculations ([ Hauschildt &| 
|Baron|2014| l may thus change the details of the spectral line profiles. Nevertheless, the cal¬ 
culations presented here demonstrate that a powerful instrument like ALMA with its high 
spatial and spectral resolution is indeed required to make full use of the diagnostically inter¬ 
esting lines at mm wavelengths. 

It also is possible that the CO pure rotational transitions could go into emission, if 
they arise in warmer material than is present deeper where the sub-millimeter continuum 
becomes optically thick. It is unlikely, however, that the sub-mm CO lines will appear in 
emission off the limb, like their thermal-IR cousins, because the sub-millimeter continuum 
forms much higher than the 5 fim counterpart, and the slant-angle effect likely will cause 
the sub-mm radiation to overwhelm the potential CO rotational transitions. In other words, 
the effective altitude where the sub-millimeter continuum would become transparent at the 
limb would be higher than the low-chromosphere region where the molecules apparently 
can survive. Nevertheless, such phenomena should be looked for, because CO has revealed 
many surprises in the past. 


2.4 Gyroresonance and gyrosynchrotron radiation 

Gyroresonance emission l |Kakinuma & Swarup|1962^|Zheleznyakov|1962a|b[ l, which is re¬ 
sponsible for the microwave emission above sunspots, is produced by thermal electrons in 
the presence of magnetic fields but is rather unlikely to be observed in the millimeter range. 
Corresponding intensity contributions would be measurable only at the lowest gyroharmon- 
ics, j = 2 - 5, which are below 30 GHz and thus outside ALMA’s frequency range even 
for the strongest magnetic fields occurring in the solar atmosphere. For instance, a mag¬ 
netic field strength of 3570 G would be required for a significant (third harmonic) emission 
contribution at a wavelength of 10 mm. 

In contrast, non-thermal gyrosynchrotron (GS) emission from energetic electrons should 
be detectable because the dominant gyroharmonics rapidly increase with electron energy. 
For example, for a relativistic electron with the Lorenz factor y, the dominant emission 
frequency is by a factor of y^ larger than the gyrofrequency, which can fall everywhere in 
the ALMA spectral domain. Such energetic particles are present in the solar atmosphere 
primarily during flares. Based on experience from the known parameter space (resolution 
and sensitivity) we therefore expect ALMA to detect non-thermal gyrosynchrotron emission 
due to high-energy electrons during times of flares or other major activity. However, with 
the new parameter space opened up by ALMA, we may also find non-thermal signatures on 
small spatial scales and short time scales in phenomena associated with intense magnetic 
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flux concentrations in network or plage regions or due to reconnection events in the quiet 
Sun, micro- or “nanoflares” (see Sect. |4.4| l. 

In general, the gyrosynchrotron absorption and emission coefficients are complicated 
functions of the energetic electron distribution, the magnetic field and the physical parame¬ 
ters of the ambient plasma ( |Ramaty|1969[ l. The standard spectrum of the gyrosynchrotron ra¬ 
diation has a bell shape, on which some harmonic structure can be seen at the low-frequency 
range (Melrose|1968 Ramaty|1969|l. The spectral peak frequency Vpeak can be specified by 
either the condition t ~ 1, demarcating the optically thick and thin regimes, or by the effect 
of high plasma density (Razin effect, Vpeak ~ 20 Although neither low- nor high- 

frequency slope is a perfect power-law, the shape of this spectrum is often characterized by 
the low- and high-frequency power-law spectral indices, aif and ahf. In the optically thick 
regime ay = 2 - 3 in a uniform source, while in the conditions of the dominating Razin 
effect, the spectrum can be much steeper, ay > 3. The source non-uniformity can make 
the observed spectrum much shallower, however. The high-frequency spectral index a/y is 
primarily determined by the fast electron energy spectrum and their pitch-angle anisotropy 
( |Melrose|1968(|Ramaty|1969[ [Fleishman & Melnikov|2003[ [Kuznetsov et al.|2011| l. There 

is no simple and accurate analytical formula reliably linking the slope of the particle distri¬ 
bution in the energy (and/or angular) domain and the radio spectral index a/y the available 
analytical approximations are too crude to be quantitatively applied. The exact numerical 
expressions ( [Melrose|1968[[^maty|1969l l are very cumbersome and computationally slow 
to be routinely used. Fortunately, fast gyrosynchrotron codes have recently been developed 
( [Fleishman & Kuznetso\j [2(^ i, which, for instance, as integral part of the 3D modeling 
tool GX Simulator ( Nita et al.|2015[ l facilitate the quick computing of the gyrosynchrotron 
emission from realistic 3D volumes (see Sect.[4.3.T]for more details). 


2.5 Polarisation and magnetic field measurements 


The chromospheric magnetic field is fundamentally important but poorly known because the 
currently available diagnostics for measuring the chromospheric magnetic field have short¬ 
comings. There are a few magnetic-sensitive lines, e.g. Call 854.2 nm and He 1 1083 nm, 
that can be used for chromospheric magnetometry but these lines are commonly charac¬ 
terized by a weak polarisation sig nal, wide formation height range and require a careful 
treatment of non-LTI0 effects (e.g., de la Cruz Rodriguez et al.|2010 1 . Most recent work on 


the weak chromospheric magnetic field has focused on the analysis of the CalllR triplet 


lines (849.8 nm. 

854.2 nm and 866.2 nm), e.g., employing the Hanle effect (e.g.. 

Carlin & 

Asensio Ramos[ 

2014[l, and on the HeI1083nm line (e.g., [Solanki et al.[[2003| 

Martinez 


Gonzalez et al.[[20l2[ [Schad et al.[[2013y Nevertheless, measuring the magnetic field is 
fundamental for advancing our understanding of the chromospheric plasma because it is 
strongly coupled to the magnetic structures. This intricate connection is vividly seen in de¬ 
tailed 3D MHD simulations (e.g., [Gudiksen et al.|201 l[|Schaffeiffieiger^^j 2006[l and has 
been widely confirmed by observations (e.g. in Ha images, De Pontieu et al.[2007 1. The role 
of the chromospheric magnetic field in dynamics, wave propagation, and (local) heating can 
hardly be overestimated. 

A method for measuring chromospheric magnetic fields by evaluating thermal brems- 
strahlung was developed by [Bogod & Gelfreikh[ ( [198^ and by [Grebinskij et ar] ( [2000[ l and 
tested on the RATAN-600 and Nobeyama Radio Heliograph (NoRH) data at short cm wave- 


Non-LTE, also abbreviated NLTE, refers to departures from local thermodynamic equilibrium (LTE). 
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lengths. The proposed technique relies on the fact that, in the general case, the longitudinal 
component Si of the magnetic field can be derived from the observed degree of circular 
polarisation V and the brightness temperature spectrum TbCv). 

The propagating of electromagnetic radiation in the magnetized solar medium is well- 
described by the magneto-ionic theory, the observable modes being the ordinary mode (o- 
mode) and the extraordinary mode (x-mode). In a magnetized plasma the absorption co¬ 
efficient and thus the opacity k for free-free emission is higher in the extraordinary mode 
than in the ordinary (/Cx > Due to the resulting dependence of the opacity on the mag¬ 
netic field the natural modes become optically thick in different atmospheric layers so that 
z(Ty_ (v) = 1) > z{To (v) = 1). Thus, at a given frequency we observe the x- and o-mode 
emission from slightly different heights with deeper (cooler) layers seen in the o-mode and 
higher (hotter) layers seen in the x-mode in the atmosphere with T increasing outwards. 

Depending on the sign of the longitudinal magnetic field component, a mode (i.e., the 
x-mode) can then correspond to either left circular polarisation (LCP) or right circular po¬ 
larisation (RCP). In the case of a magnetized atmosphere with a temperature gradient, the 
temperature difference between these layers produces then an observable net circular po¬ 
larisation. The stronger the magnetic field is, the larger is the difference of the formation 
heights of the x- and o-modes, resulting in a proportionally stronger net polarisation. See 
Fig. and the corresponding text in Sect. |4.3.1| for simulated examples, illustrating how 
the degree of net circular polarisation is determined from brightness temperature maps for 
left and right circular polarisation. 

|Bogod & Gelfreikh] ( |1980[ l and [Grebinskij et al.| l |2000[ l demonstrated that the slope of 
the (observed) free-free emission brightness temperature spectrum (T^(v)) 

_ d\og{Tb(v)) 

^ d\og(v) ^ 

carries information about the gas temperature gradient in the radiation forming layers in the 
atmosphere and controls the degree of the circular polarisation P, such that 

V = ^— cose , (13) 

V 

where v is the observing frequency, vb is the (electron) gyrofrequency 

1 eB 

vb=- - =2.8xlO®HzxB , (14) 

In meC 

and 9 is the angle between the line of sight and the magnetic field vector. In the simple 
case of an optically thin homogeneous atmospheric slab with a constant gas temperature 
the quadratic dependence of the observed brightness temperature on the observing fre¬ 
quency V (cf. Eq. [^) results in a local slope f = 2. The circular polarisation, P, is then 
simply (cf. |Kundu[l965[ l 

r 2— cose . (15) 

V 

This expression is valid for quasi-longitudinal propagation, which prevails unless 9 is close 
to 90 degrees. Putting numerical values we get: 

P (%) 1.85 X 10“^T (mm) 5 (G) cos0 (16) 

which gives P - 2% at 1 mm for a field of 1 kG as it occurs above sunspots (see, e.g., 
|Solanki|2003[ and references therein). In the optically thin case the longitudinal component 
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of the magnetic field S/ can thus be found from observations at a single wavelength, but 
spectral observations are still desirable to confirm the free-free nature of the emission and 
its thickness. Detailed modeling for the quiet Sun (see Sect. |4.2.^ predicts for the optically 
thick chromosphere a circular polarisation of 0.05% at 1 mm and thus gives an estimate of 
the required accuracy for polarisation measurements, although the polarization can reach a 
few % above sunspots (see Sect. |4.3.T] l. 

The expression in Eq. offers an elegant means of measuring the longitudinal com¬ 
ponent of the magnetic field using the final expression (together with Eq. ( |14[ l) 

B, = |B|cos0= ^ (2.8 X 10®)“* Hz“' (G) . (17) 

The restored magnetic field refers to the height where the emission at a given frequency 
is formed so that the magnetic field at different heights can be measured by changing the 
frequency (see Sect. |3.3.4| for details for the case of ALMA). 

It should be noted that the approach outlined above implicitly assumes a constant mag¬ 
netic field within the layer between Tx = 1 and Tq = 1 and thus neglects the complex mag¬ 
netic field structure on small spatial scales as it is known from contemporary high-resolution 
observations, although, depending on the local magnetic field strength, the height difference 
between Tx = 1 and Tq = 1 can be small enough to allow for such an assumption. Please 
refer to Sect. |3.3.4| for a description of ALMA’s potential for measuring the chromospheric 
magnetic field based on the method introduced here. 


3 Observations of the Sun at millimeter and submillimeter wavelengths 


3.1 Solar radio observations in the pre-ALMA era 


In spite of the important progress during the last forty years, the mm and sub-mm range 


has remained by large unexplored, except for single dish observations (e.g. Trottet et al.| 
201 1[ and references therein) and some observations with BIMA l|Loukitcheva et al.|2006[ 


White^et al J2006^ . An important step forward after the first interferometric observations 


of Kundu (1959 1 was the use of non-solar aperture synthesis instruments for solar observa¬ 
tions. The first high resolution radio imaging observations of the Sun were obtained with 
the Westerbork Synthesis Radio Telescope (WSRT) about forty years ago, first of active re¬ 
gions at 6 and 20 cm ( Kundu & Alissandrakis|1975] |Chiuderi Drago et al.|1977 and bursts 


( [Alissandrakis & Kundu|1978| l and subsequently of the quiet Sun (|Kundu et al. 1979| l. De¬ 

tailed modeling of the sunspot associated emission ( [Alissandrakis et al.|1980 1 gave impor¬ 
tant information both on the emission mechanism and the structure of the atmosphere, in 
particular the magnetic field structure in the transition region and low corona, while model¬ 
ing of the microwave emission from flaring loops ( [Alissandrakis & Preka-Papadema 1984^ 
made possible the association of their observable characteristics with the properties of the 
trapped energetic electrons. Quiet Sun observations showed a marked similarity to the chro¬ 
mospheric network as it was known from observations in other wavelength ranges. Ever 
since, solar observations with the Very Large Array (VLA) have provided increased spatial 
and temporal resolution and extended the frequency range up to 15 GHz (2 cm, see, e.g., 
[Bastian et al.|1998t|Lee|2007[|Shibasaki et al.|2011t[Gary & Keller|2004| l. In parallel with 
the occasional very high spatial resolution observations with the WSRT and the VLA, so¬ 
lar two-dimensional images with moderate resolution were taken daily with the Nobeyama 
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Radioheliograph (|Nakajima et al.||1994 

at 17 and 34 GHz (1.76 and 0.88 cm), the Met- 

sahovi Radio Observatory dVrsnak et al. 

1992[l at 37 GHz (8 mm), the Siberian Solar Radio 

Telescope (SSRT, Grechnev et al. 2003a i at 5.74 GHz (5.2 cm), the Nan 9 ay Radioheliograph 

(|Kerdraon & Delouis|1997 

1 at frequencies between 150 to 450 MHz (67 to 200 cm), together 


with the one-dimensional, high frequency resolution data of RATAN-600 l[Bogod| 2011[|Bo-| 
|god et al.|2011[ ) from 3 to 18.2 GHz (1.65 to 10 cm). The Owens Valley Solar Observatory 
(OVSO) has also provided valuable data. 


3.2 The Atacama Large Millimeter/submillimeter Array 

ALMA is an international partnership between Europe, North America, and East Asia in 
cooperation with the Republic of Chile and minor partners to build and operate a millime¬ 
tre/submillimeter interferometer on the Chajnantor plateau in the Chilean Andes at an al¬ 
titude of 5000 m. The 66 antennas of this aperture synthesis telescope are arranged in two 
major groups: The 12-m Army and the Atacama Compact Array (ACA). The 12-m array 
consists of 50 movable antennas with 12 m diameter. The ACA, also known as “Morita Ar¬ 
ray”, combines 12 antennas with 7 m diameter for interferometry (the ACA 7-m Array) and 
4 antennas with 12 m diameter for single dish observations surrounding them (also referred 
to as the ACA Total Power (TP) Array). The ACA antennas are installed in a very com¬ 
pact configuration, while the more extended 12-m Array can be rearranged to form compact 
or more widely spread configurations with baselines, i.e. distances between the individual 
antennas, of up to 16 km. The ACA TP Array provides information equivalent to an in¬ 
terferometer with baselines from Om up to 12 m. The ACA 7-m Array samples baselines 
from 8.5 m to 33 m, bridging the baseline sampling gap between the 12-m Array and the 
ACA TP Array. 

The maximum spatial resolution is determined by the observing frequency and the 
longest baseline of the 12-m Array, and can reach a few milliarcseconds. Images with high 
fidelity at the maximum resolution, however, can only be expected if the brightness dis¬ 
tribution is dominated by bright compact emission. For complex images of chromospheric 
structures, where the detected emission fills the whole beam, the effective spatial resolution 
may be limited to a few 100 milliarcsec at a wavelength of 1 mm. This will be quantified 
more accurately on the basis of further observations (see also Sect. |3.3.T] l. 

The ACA TP array is an important part of the overall array because without it ALMA 
would suffer from a problem common to all synthesis instruments: The lack of short base¬ 
lines and consequent short spatial frequencies would make extended sources invisible. The 
maximum recoverable size, as defined in the ALMA technical handboolQ is 

0.64 , 37100 

omax = -rad = -arcsec (18) 

^inin ^ 

where L,„i„ is the minimum baseline in meters and v the observing frequency in GHz. For 
l^min = 8.4 m (7-m array, NS), this gives &,„ax ~ 44 arcsec at 100 GHz. This simply means 
that, without the ACA TP array, the quiet sun background and also large sunspots would 
be invisible although still fine strucfures in the chromospheric network and in solar flares 
would be detectable. See Sect. 13. 3. Tl for a discussion of interferometric observations of the 
Sun and the role of the array components. 


7 


http://almascience.eso.org/proposing/call-for-proposals/technical-handbook 
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The instantaneous field of view (FOV) is conventionally given by the full-width-half- 
maximum (FWHM) of the main lobe of a single antenna and depends on the wavelength A, 
antenna diameter D and its illumination function. The actual ALMA 12-m antennas have a 
measured primary beam FWHM, which scales with wavelength approximately as 

6»~ 1.13 X - ~ 19" X for D= 12m. (19) 

D 1 mm 

Therefore ALMA can observe only a small fraction of the Sun at a given time. However, 
the total FOV can be increased by rapid mosaicing (see Fig.|20|> or by “on-the-fly” interfer¬ 
ometric scanning, i.e. by continuously changing the pointing. 

Every antenna will eventually be equipped with receivers covering ten frequency bands 
in the range from 35 GHz to 950 GHz, corresponding to wavelengths from 8.6 mm to 0.3 mm. 
ALMA can observe with one receiver at a time but the hardware is designed to switch be¬ 
tween up to three pre-selected receivers in about 1.5 s. The exact achievable band-to-band 
cadence is determined by various overheads and has yet to be determined exactly. 

Each of the ALMA bands has two polarisation channels with an instantaneous total 
bandwidth of 8 GHz per polarisation. These are down-converted into 2 GHz wide base¬ 
bands, digitized and sent to a Correlator. The Correlator for the 12-m array can operate 
in two main modes, TDM (Time Division Mode) for low resolutions wideband continuum 
observations, and FDM (Frequency Division Mode) for higher spectral resolutions. In FDM 
mode, each 2 GHz baseband is subdivided by digital filtering by up to 32 sub-bands, each 
with 62.5 MHz bandwidth (usable bandwidth 58.6 MHz to avoid edge and aliasing effects). 
These sub-bands can be individually tuned across the baseband and allow to control the 
analysed bandwidth and thus the frequency resolution. The maximum number of frequency 
channels in FDM dual polarisation mode is 3840. Thus, if all sub-bands are stitched to¬ 
gether, the full baseband will be covered with spectral channel spacing of 488 kHz. On the 
other side, if all Correlator resources are used to analyse one sub-band for ultimate spectral 
resolution, the corresponding channel spacing is 15.3 kHz. The absolute minimum channel 
spacing achievable is 3.8 kHz with the 29.3 MHz bandwidth and in single polarisation mode. 
The resulting spectral resolution is 7.6 kHz (larger than the channel spacing), which corre¬ 
sponds to a velocity resolution of 0.02km s“* at T = 2.73 mm (v = 110 GHz). Averaging of 
spectral channels is available in multiples of 2, to reduce data rate in FDM mode. The TDM 
mode provides up to 128 channels per baseband (~ 120 usable channels due to truncation of 
the edge channels in offline data processing), covers the full baseband bandwidth and pro¬ 
vides channel spacing of 15.6MHz in dual polarisation mode. This mode has an advantage 
of lower data rate and faster dump time, so it is used for all observations where high spectral 
resolution is not required. For more information, see the ALMA Technical Handbook^. 

While ALMA is a general purpose telescope, provisions have been made to support so¬ 
lar observing. First, the surface of the antennas has been roughened to ensure that it scatters 
much of the flux at visible/IR wavelengths so that the receivers are not damagecQ Further¬ 
more, the detectors of ALMA are very sensitive, while the solar radio flux is very strong. 
At the same time most of the instrumental calibration steps, like pointing, focus corrections, 
and phase calibrations have to be done using much weaker radio sources, e.g. quasars. When 
pointing at the Sun, it is necessary to reduce the high flux at millimeter wavelengths in order 
to avoid saturation of the receivers. This is done by either inserting solar Alters in the optical 
path, or by electronic detuning of the detectors. Test observations have shown that obser¬ 
vations of the quiet Sun and moderate solar flares can proceed without the use of the solar 


Please refer to the ALMA websites for more technical information (http://www.almaobservatory.org/ 1 . 
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attenuators, via the detuning alternative method (although usage of the filters remains an op¬ 
tion). This may also be possible for observing major flares given that the stronger flux in big 
events comes from a bigger area; thus, even the most energetic flares will still look ’moder¬ 
ate’ within the restricted ALMA FOV. The exact procedure has yet to be tested on stronger 
flares. The receiver sensitivity is degraded in the detuned mode, with the noise tempera¬ 
ture typically in the range 100-1000 K, depending on the required attenuation. However, 
on the solar disk, the sensitivity will be still limited mostly by temporal and spatial fluc¬ 
tuations (e.g., p-mode oscillations and propagating waves) of solar origin - “signal” rather 
than “noise” in some sense. The test data from late 2014 have convincingly confirmed this 
conclusion. 


3.3 alma’s diagnostic potential for the Sun 

The technical specifications summarised in Sect. |3.2| promise impressive diagnostic oppor¬ 
tunities for solar observations with ALMA. In the following, we discuss the advantages and 
challenges of observing the Sun with an interferometer and highlight key capabilities. 


3.3.1 High-resolution interferometric imaging at millimeter wavelengths 

The Sun is a very difficult object to observe with an interferometer, and it can be argued 
that the solar chromosphere is the most difficult target of all. The reason is that the Sun 
fills any field of view with emission on a wide range of spatial scales, from the large-scale 
background down to the smallest scales of magnetic field structures, but interferometers are 
designed to sample a fixed range of spatial scales corresponding to the range of baseline 
lengths in the arra)|^ As an example, a single observation with the Very Large Array (VLA) 
covers about a factor of 20 in spatial scales. Synthesis of a large aperture by this technique 
requires reasonably dense sampling of spatial scales for successful calibration and robust 
imaging, so one cannot simply impose an arbitrary range of spatial scales and achieve a 
successful observation ( |Taylor et al.||199^ . Arrays such as the VLA and ALMA can get 
around this limitation of a single configuration by observing the same object with multiple 
configurations, each covering complementary ranges of spatial scales: thus the VLA can 
sample a factor of order 750 in spatial scales by observing the same object in all 4 available 
configurations. However, since configuration changes for the VLA and ALMA typically 
take at least a week, this technique for expanding the range of spatial scales only works for 
objects that do not vary on corresponding timescales, and thus is completely useless for the 
Sun. 

We are therefore limited to the range of spatial scales present in a single ALMA configu¬ 
ration^. When only the 12 m antennas are used, the minimum spacing in the interferometer 
array is around 14 m and this limits the largest spatial scale sampled in all configurations. 
Most of the ALMA configurations currently in use have a range of about 20 in spatial scale, 
but this can be enhanced by including the array of 7 m antennas for which the minimum 
spacing is around 8 m. As mentioned earlier, for bright compact sources such as solar flares 
where emission is dominated by sources with a limited range of spatial scales, this range of 
spatial scales is quite adequate for excellent imaging (see Sect. |4. 4. T] regarding imaging of 

® The angular scale on the sky sampled by an interferometer formed from a pair of antennas is proportional 
to the wavelength divided by the physical separation of the antennas projected onto the plane perpendicular 
to the line of sight to the source. 
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flares). This is also true of VLA observations of coronal sources, which are typically at least 
several times brighter than the background solar emission. However, the solar chromosphere 
at millimeter wavelengths is composed of emission on all spatial scales, with contrasts of 
up to lOOOK brightness temperature that must be imaged against a continuously changing 
background of order 4000 - 7000 K. Thus the spatial structure to be imaged by the interfer¬ 
ometer is much weaker than the background. Since it has components on all spatial scales 
from subarcsecond out to the field of view of the dishes, imaging is challenging. 

The issues involved in imaging the solar chromosphere at millimeter wavelengths were 
discussed by [White et al.| ( (2006[ l in the context of the BIMA observations shown in Fig.[^ 
The critical issue is sampling the largest spatial scales of the emission: if this is not sam¬ 
pled adequately then images can be dominated by spurious small-scale features. [White et 3T] 
( |2006| l addressed this problem by using the most compact configuration available at BIMA 
and by deconvolving using the maximum entropy method (MEM) to take advantage of the 
information on large spatial scales present in the primary beam response of the dishes (| Ekers| 
|& Rots|1979| l. However, since short baselines are essential to sample the large-scale struc¬ 
ture and interferometers usually sample a limited range of spatial scales, one usually has to 
sacrifice spatial resolution to achieve successful imaging. 

ALMA possesses three critical features that should allow successful imaging of the 
chromosphere with much better spatial resolution than previously, (i) The larger number 
of dishes (up to 62 for interferometry) provides much more information in each integration. 

(ii) The second advantage is the availability of the AC A array of 7 m dishes: as long as it is 
in a compact configuration it can provide the shorter spacings lacking from the wider 12-m 
array that provides high spatial resolution. Just how wide a range of spatial scales can be 
achieved in this fashion remains to be demonstrated: it may be that we will not be able to 
successfully image the Sun at the highest spatial resolution of ALMA because the amount of 
information needed to reconstruct a filled field of view (up to 10^ resolution elements) is not 
adequately sampled. The fact that the field of view of the 7 m dishes is almost twice as large 
as the field of view of the 12 m dishes may also cause problems in deconvolution. However, 
preliminary experiments using CARMA observations with both 6 m and 10 m dishes indi¬ 
cate that successful imaging is still possible (White, Loukitcheva & Solanki, in preparation). 

(iii) The other feature of ALMA that will be a major factor is the availability of the 12 m 
TP dishes. These dishes have been designed with the ability for scanning the sky at very 
fast rates, and have already been shown to be capable of imaging the full solar disk in 8 min 
(Hills & Phillips, private communication). This means that the smaller targeted interferom¬ 
eter regions can be scanned, and the zero spacing flux recovered, on timescales shorter than 
1 min, fast enough to track coherent 3-minute oscillations in the chromosphere. The images 
from the TP measurements in combination with the interferometer data, with the 7-m array 
providing intermediate spatial scales, is expected to improve the imaging ability of ALMA 
dramatically. In addition the single-dish data will provide absolute measurements of the to¬ 
tal solar brightness temperature, including the background disk level that is filtered out by 
the interferometer, and thus provide interferometer maps with the correct temperature scal¬ 
ing. Fast-scanning of the whole solar disk (or later a small FOV) with the TP antennas will 
most likely be offered as an integral part of interferometric observations, aiding the image 
reconstruction and providing context information. 

As briefly discussed before, the effective spatial resolution might be much lower when 
reconstructing an image for an extended area source like the Sun from a limited number of 
baselines compared to what is theoretically possible with very long baselines for separating 
two point sources. The effective spatial resolution of a reconstructed image is determined by 
the coverage in the spatial Fourier space (the u-v-space) and the corresponding synthesised 
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point spread function of the combined array. The spatial resolution depends then on the ar¬ 
ray configuration, the observed wavelength and the position on the sky (due to the elevation 
angle of the antennas and the resulting foreshortening of the effective baseline). Current esti¬ 
mates based on numerical simulations of ALMA’s imaging capabilities (e.g. |Wedemeyer &| 
|Parmer|2015| l imply a spatial resolution of 0”.3 x T/1 mm to 0”.4 x T/1 mm although it has 
to be quantified more accurately on the basis of actual observations. Reconstructed images at 
the shortest wavelengths, i.e. 1 mm and below, will have an effective spatial resolution close 
to the best possible with today’s optical solar telescopes. Due to the increase of continuum 
formation height with wavelength (see Fig.[^, ALMA samples the low chromosphere with 
higher spatial resolution than the upper chromosphere. However, it is important to empha¬ 
sise that each observation still contains information on even smaller scales, produced by the 
longest baselines. There will certainly be scientific applications for which the visibilities, 
i.e. the measured signal on a baseline, can be used directly without the need of image recon¬ 
struction, thus utilising information on much smaller spatial scales. In principle, this is true 
for all studies that are carried out in spatial Fourier space. A potential application could be 
the determination of the spectral characteristics and levels of MHD turbulence in the solar 
chromosphere, which is an important question in solar physics. Nevertheless, the anticipated 
spatial resolution of reconstructed ALMA maps will be substantially higher than everything 
that has been possible at millimeter wavelengths before (cf. Fig.[^. 

3.3.2 High cadence 

In interferometry, the integration time is determined by two factors: the sensitivity and the 
coverage of the u-v plane. Sensitivity is not a problem because the Sun is a very bright 
source at sub-millimeter and millimeter wavelengths and because ALMA is a very sensitive 
instrument. For interferometric observations of simple structures with strong contrast, such 
as bursts, snapshot mode (i.e., instantaneous measurements) and thus ultra-high cadence 
imaging should be possible with ALMA. The situation is more complicated if the observa¬ 
tional target has a complicated structure, which fills the primary beam. As mentioned above, 
in snapshot mode, ALMA samples the u-v plane at discrete points. If the sampling is uni¬ 
form, in steps of Ab wavelengths, the beam pattern is repeated at intervals of XjAb radians. 
Even if the sampling of the u-v plane is not uniform, as is the case with ALMA, confusion 
will still arise if the field of view has a complex structure as, e.g., in the quiet sun. The 
u-v coverage required for resolving the structure puts some constraints on the achievable 
cadence, usually leading to a compromise between spatial and temporal resolution, which 
has to be adjusted according to the requirements of individual scientific applications. It is 
possible to take advantage of Earth’s rotation for increasing the number of sampled points in 
the u-v plane (then in small arcs), which then requires longer integration times. Integrating 
for several minutes or even up to 40 min, as it has been done during the fifth Commissioning 
and Science Verification campaign in December 2014, is reasonable for studies of slowly 
evolving features but is prohibitive for observing fast evolving phenomena. For the latter, a 
time resolution of less than 1 min is required. 

A cadence well below 1 s should be achievable for snapshot imaging, which enables the 
study of high-frequency waves in the solar chromosphere. These waves are interesting re¬ 
garding their role in the energy transport and heating (see Sects. [43] and |4.6[ l but also provide 
a means of probing the atmospheric structure (“chromospheric seismology”). For instance, 
eclipse observations were previously used to detect waves with periods of 6 s in low coronal 
loops ( [Williams et al.|2002| l. Such coronal loops will be observable also with ALMA, but 
much more often and for longer time sequences than possible during solar eclipses. It is also 
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expected that the miniature loops found by [Peter et SI] ( |2013| |2014[ l have high frequency 
oscillations, because of the smaller length scales (but wave speeds comparable to regular¬ 
sized active region loops). ALMA will provide the temporal and spatial resolution needed to 
observe such oscillations. Moreover, the high cadence ALMA observations will be useful in 
observing oscillations observed in flare light curves, the so called quasi-periodic pulsations 
(QPP, for a review see e.g. [Nakariakov & Melniko^|2009| l. Please refer to Sect. [4.4.^ for 
quasi-periodic pulsations and Sect. |4.6| for more details on chromospheric waves. 

Apart from high-frequency waves, the relevant temporal scales for phenomena in the 
solar chromosphere are typically on the order of 10 s and in some cases even longer. Multi¬ 
band (and thus multi-layer) data cubes could be obtained when switching receiver bands 
quickly with two to three of them being ready for instant use. It is not entirely clear yet what 
effective band-to-band cadence can be achieved and if it would be fast enough for resolving 
the shortest timescales. 

The achievable cadence for mosaicing depends on many factors such as the required 
FOV and thus the number of pointings, the required integration time per pointing, and the 
required spectral resolution, which determines the correlator mode and thus its minimum 
dump time (with the FDM mode being generally slower, see Sect. |3.2| i, and other over¬ 
heads like the antenna slew rate. On-the-fl)!*^ interferometry (OTFI) is probably faster for 
low time per pointing but has yet to be commissioned. The simultaneous use of sub-arrays 
acting in different frequency bands is also possible in principle, but not in the initial solar 
observations. The exact setup of mosaicing observations and the resulting cadence can and 
should be adjusted to the requirements of the considered scientific aim, finding an optimum 
compromise between FOV, temporal and spectral resolution. 

3.3.3 Spectral capabilities 

The high achievable spectral resolution and the large number of simultaneous frequency 
channels enables ALMA to record precise spectra. Recombination lines (i.e. Rydberg tran¬ 
sitions, see Sect. |2.2| l and molecular pure rotational transitions (Sect. |2.3| l could thus be 
utilized as powerful complementary diagnostic tools with ALMA. Among the ionic species 
with lines in the wavelength range accessible to ALMA are, for instance, H I, C III, OIV, 
O V, O VI, Ne VII, Si XI , Fe XV. While the continuum radiation is mostly formed in the 
chromosphere, some of the recombination lines originate in the corona, thus adding coronal 
diagnostics to ALMA’s tool kit. Some of these recombination lines probably are severely 
broadened, which has made detection with previous instruments difhcult. The high spec¬ 
tral resolution and high sensitivity of ALMA could finally lead to their detection. Carbon 
monoxide is probably the most promising molecule with its pure rotational spectrum in 
the ALMA wavelength range, especially because it provides a local kinematic and thermal 
probe for the dynamic solar chromosphere. Observations of CO also will benefit from the 
high spatial and temporal resolution because the CO number densities are affected by hot, 
propagating shock waves on small spatial scales. 

3.3.4 Polarization capabilities and chromospheric magnetic field measurements 

ALMA’s capability of measuring the polarisation state at multiple wavelengths will provide 
a powerful diagnostic tool for the magnetic field over a subsfantial range of chromospheric 

In contrast to a sequence of discrete pointings, the antennas keep moving during an “on-the-fly” obser¬ 
vation, thus changing the pointing continuously and scanning the target along predefined trajectories. 
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heights. For a magnetized atmosphere with a temperature gradient, like it is the case for 
the Sun, the temperature difference between two layers produce an observable net circular 
polarisation (see Sect. [ 23 , which can be measured by ALMA. As mentioned in Sect. |2.5[ 
the height at which the magnetic field is determined varies with frequency. Strictly speaking, 
the difference in formation heights between the o- and x-modes should be taken into account 
but models show that the difference is small (see Sect. |4.2.^ . The contribution function of 
the continuum intensity displayed in Fig. can thus also be understood as the effective 
height where the magnetic field is probed. In the case of ALMA, it gives access to magnetic 
field information from the temperature minimum layer up to the transition region over the 
ALMA spectral range. 

The technical requirement for the accuracy of the polarized flux^’* is 0.1% of the to¬ 
tal intensity, after calibration. With the current on-axis calibration scheme, it is feasible to 
achieve this accuracy for a compact source at the beam centre. Polarisation imaging ac¬ 
curacy degrades off-axis, up to 0.3% for angular source sizes of FWHM/3 of the primary 
beam. To achieve 0.1% accuracy for extended sources, the off-axis instrumental polarization 
calibration has to be performed. As a result, high-precision multi-frequency observations of 
the polarized emission with ALMA should provide a wealth of information for the restora¬ 
tion of the longitudinal component of the magnetic field over a wide range of heights in the 
chromosphere. More detailed studies based on state-of-the-art numerical models and high- 
resolution observations (including ALMA) are currently carried out in order to investigate 
the limits and reliability of the approach. Based on these results, magnetic field measure¬ 
ments in active regions (see more details in Sect. |4.3.T] i will clearly be possible while such 
measurements for quiet Sun regions (more details in Sect. |4.2.3| l are at least very likely to be 
feasible given ALMA’s high sensitivity. Polarization measurements will also be extremely 
valuable during flares, which will put strong constraints on the pitch-angle anisotropy of 
accelerated electrons and will be able to distinguish electron and positron contributions (see 
Sect. 14.411. 


3.3.5 ALMA observations in co-ordination with other observatories 

To make the most of ALMA’s unprecedented spatial resolution, and also to place its inten¬ 
sity and polarisation measurements into a broader context, observations at other wavelengths 
are essential. In particular, several ground-based optical observatories with l-t-meter-class 
apertures, such as Big Bear Solar Observatory’s L6-m New Solar Telescope (NST, [Goode &| 
|Cao|2012) and the European GREGOR L5-m telescope ( [Schmidt et al.|201^ , are producing 
diffraction-limited observations in the visible and near-infrared. In addition, other ground- 
based observatories can provide important context data, such as the Swedish 1-m Solar Tele¬ 
scope, in particular known for its chromospheric observations with the CRISIp]instrument, 
Themis, VTT and Chrotel on the Canary Islands, Kanzelhohe Solar Observatory in Austria, 
Hvar Observatory in Croatia, and also the Daniel K. Inouye Solar Telescope (DKIST, for¬ 
merly the Advanced Technology Solar Telescope, ATST, |Keil et al.|2011| currently under 
construction, possibly operational by 2019/2020) and the future European Solar Telescope 
(EST, [Collados et al.|201^ . Eor instance, the DL-NIRSP (Diffraction limited Near Infra- 
Red Spectropolarimeter) instrument on DKIST will provide valuable chromospheric (and 
possibly coronal) magnetic field measurements with a spatial resolution of ~ 60 km in the 
NIR and down to ~ 25 km for Ho' (in the most ideal case). The New Solar Telescope pro- 

* * The Swedish 1 -m Solar Telescope will soon he equipped with the CHROMIS instrument which will have 
a narrow-band Call K filter suitable for chromospheric observations. 
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Fig. 7 TiO band (705.7 nm) continuum image from NST of a sunspot and nearby granulation on 2014 Jul 31 
(courtesy of Vasyl Yurchyshyn, BBSO). The image has been resampled with a pixel size of 0.07 arcsec/pixel, 
whereas the original image from the camera had 0.”034 pixels. The angular resolution is close to the diffrac¬ 
tion limit of the telescope, which is 0.”09. The circle illustrates ALMA’s primary beam at wavelengths around 
1 mm. 


duces images, which have an angular resolution close to the diffraction limit of the telescope, 
i.e. 0.”09 at 705.7 nm (in the TiO band, see Fig.|^ and 0.”16 at 1560nm. With the ability to 
measure magnetic helds on this sub-arcsecond scale with the NIRIS instrument ( |Cao et ah] 
|2012[ l and others, these observations of both the deeper (at 1,560 nm) and higher layers (e.g., 
in Ha or He 11083 nm) will provide an ideal complement to the ALMA high-resolution in¬ 
terferometric observations of the quiet Sun and active regions. At Big Bear Solar Observa¬ 
tory, another complementary instrument to ALMA is Cyra, a cryogenic Near Infrared (NIR) 
echelle spectrograph ( [Goode & CaopoO] ! capable of spatially-scanned, high spectral reso¬ 
lution slit imaging out to 5000 nm. This includes the CO molecular lines at 4,667 nm, which 
can directly support ALMA’s exploration of the coolest parts of the solar atmosphere (see 
Sects. [4~^ and |2.2] l. 

At GREGOR, the GRIS (GRegor Infrared Spectropraph) instrument ( [Collados et al.j 
[2012[ l is providing measurements close to the diffraction limit of the velocity and magnetic 
field in the upper chromosphere in the He 11083 nm line, as well as in the photosphere at 
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Fig. 8 Synthetic maps based on a 3D numerical simulation of an enhanced network region with magnetic 
loops connecting patches of opposite polarity: (a) UV continuum close of the Mg II k line showing the granu¬ 
lation at the bottom of the photosphere, (b) Vertical magnetic field component at the bottom of the photo¬ 
sphere. (c) Brightness temperature at 1 mm. (d) Brightness temperature at 3 mm. (e) Intensity in the blue line 
wing of the Mg II k line, (f) Intensity in the line core of the Mg II k line. The contrast is written in the upper 
right comer of each panel with exception of panel b for which the standai'd deviation is specified. The circle 
marks the primary beam of ALMA at T = 1 mm but is not shown for the 3 mm map because, in that case, 
the primary beam is three times in diameter and is thus larger than the depicted region. The numerical model 
jCarlsson et al.|2Q16) was computed with the Bifrost code jGudiksen et al.|2011| see Sect. |4.2.'^ for more 
details). The intensity images for panels a,e and f are produced with the radiative transfer code MULTI3D 
by|Leenaarts et alT] i |2013^ , while the millimeter maps in panels c-d were calculated by |Loukitcheva et aL\ 
j2Q15a^. The magnetic field component shown in panel b is taken directly from the numerical model. 
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1560 nm, thus nicely straddling the layers. The ability to measure the full magnetic vector 
by GRIS will complement the longitudinal field measurements made by ALMA and will 
greatly facilitate the interpretation of the ALMA magnetic field data. Conversely, the very 
high cadence achievable with ALMA with greatly extend the capabilities of sampling the 
dynamics of the chromospheric magnetic field. Also, the Gregor Fabry Perot Interferometer 
(GFPI, |Puschmann et al.|20T2) l will provide high-resolution images and magnetic field data 
in various photospheric and chromospheric spectral lines. 


Simultaneous campaigns with other ground-based observatories would be limited by 
their geographic locations and resulting time differences with respect to the ALMA site, 
although some scientific applications may not require exact simultaneous observations. In 
contrast, space-borne telescopes can provide simultaneous data around the clock (except 
for eclipse periods) and offer complementary diagnostics, which are not accessible from 
the ground (e.g., EUV). The most relevant space missions in this respect are the Solar Dy- 
namic s Observatory (SDO,|Lemen et al.|2012 1 , the Interface Region Imaging Spectrograph 
(IRIS, [De Pontieu et al.|2014b i, Hinode (Kosugi et al.|2007|l, the Reuven Ramaty Fligh En¬ 


ergy Solar Spectroscopic Imager (RFIESSI, 


Observatory (SoHO,pomingo et al.|1995a 


Lin et al.|2002 1 , and the Solar and Heliospheric 
b[ [Harrison et al.|1995[|Kohl et al.|1995l l. Most 


instruments have a limited field of view just like ALMA itself and thus require detailed plan¬ 
ning prior to the observation. SDO, on the other hand, observes the whole Sun at all times 
and offers amongst many things vector magnetograms of the photosphere recorded with the 
Helioseismic and Magnetic Imager (HMI, jScherrer et al.|20r2l l, as well as filtergrams for 
diagnostics isolating many layers in the solar atmosphere (e.g., with the Atmospheric Image 
Assembly, AIA, |Lemen et al.|2012| l. Furthermore, for some applications, coordination with 
other radio telescopes or arrays could be desirable (see Sect. |4.3.T] for a short list of relevant 
radio facilities). It should also be mentioned that the Spectrometer Telescope for Imaging 
X-rays (STIX, [Benz et al.|2012b] l on board Solar Orbiter will be the only high-resolution 
X-ray instrument beyond the lifetime of RHESSI. 

For instance, the combination of ALMA and IRIS promises a set of complementary 
plasma diagnostics with great scientific potential. Since its launch in June 2013 as a NASA 
Small Explorer, IRIS has already provided fascinating new insights into the dynamics of the 
chromosphere and possible new heating mechanisms (e.g., De Pontieu et al.|[2014a[[Tian| 


|et al.|2014a[[Testa et al.|2014a[[Peter et al.|2014[|Hansteen et al.|2014 1 . It has also become 

clear, however, that the interpretation of the spectra recorded in the optically thick Mg II h 
and k lines is very difficult and will require sophisticated new analysis tools in combination 
with state-of-the-art forward modeling. An example for a complementary ALMA-IRIS data 
set is illustrated based on a numerical model snapshot in Fig. It combines (synthetic) 
ALMA maps at different wavelengths, thus probing different layers of the chromosphere, 
with important context data (here the intensity and magnetic field in the low photosphere) 
and complementary other diagnostics, namely (synthetic) IRIS maps in the Mg II k line at 
a wavelength of 279.6nm. Beside the differences in formation height, the Mg Ilk line core 
intensity is a highly non-linear measure of the atmospheric conditions. Both the ALMA 
3 mm map and the IRIS Mg II k line core map show an imprint of the magnetic loops in the 
numerical model but the contrast is dTb/ < Tb >= 21.9% for the 3 mm map and dlj < 
I >= 179.2 % for the Mg Ilk line core map. Much of the contrast in the Mg map is due to 
very localised brightenings, which can outshine the surrounding region by a factor ~ 10. 
This extreme effects make a physical interpretation difficult, especially when compared to 
the linear behaviour of the ALMA maps, but nevertheless probes different aspects of the 
chromospheric plasma conditions and thus contributes to comprehensive diagnostic data 
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set. While ALMA alone will without any douht produce excellent data, its scientific output 
can be significantly enhanced in combination with other observatories. 


4 Potential solar science cases for ALMA 

This section outlines some potential science cases for ALMA. In this respect, numerical 
simulations of the solar chromosphere can play an important role for the planning, optimiza¬ 
tion and interpretation of observations with ALMA. Synthetic brightness temperature maps 
calculated from numerical models can be used to simulate what ALMA might observe. Dif¬ 
ferent instrumental configurations can be tested and adjusted to the scientific requirements, 
finding the optimal set-up for different scientific applications. In addition, coordinated ob¬ 
servations of other chromospheric diagnostics, say at optical, infrared (IR) or ultraviolet 
(UV) wavelengths, can address open questions whose solution benefits from a broad-based 
approach to which ALMA could contribute significantly. 


4.1 Central questions in solar physics to be addressed with ALMA 

ALMA will be able to advance our knowledge of the chromosphere in all its different 
flavours ranging from quiet Sun regions to flares. In particular, substantial progress can 
be expected regarding the following Grand-Challenge topics in contemporary solar physics. 


Coronal and chromospheric heating is a long-standing problem in modern astrophysics 
(see, e.g., |Ulmschneider et al.|19^ . Despite much progress during the past decades, several 
fundamental issues concerning the dynamics of the solar chromosphere are still unresolved. 
For instance, the verdict as to what mechanisms contribute most to heating the network 
and the internetwork chromosphere is still open. Equally, the source of the so-called “basal 


flux” (|Schrijver|1987[l is still debated. Is it due to acoustic waves or processes connected 


to the Sun’s magnetic field or both (see, e.g., Schroder et al. 

2012||Cuntz et al.|2007HJudge| 

|et al.|2003[ [Judge & Carpenter|1998[[Buchholz et al.|1998 

Schrijver|1995 and references 


therein)? Or in other words, are those heating processes of the alternating current (AC) or the 


direct current (DC) type? Magnetohydrodynamic waves (e.g., |De Pontieu et al.|2007| l can 
in general contribute directly or indirectly to heating the atmospheric plasma, i.e. through 
perturbation of the magnetic field, resulting in damping of the waves and the associated 
release of magnetic energy. Do the problems require a fully kinetic model approach? And 
is there a magnetically driven basal flux connected to the presence of an ubiquitous, local, 
turbulent dynamo ( |V6gler & Schussler|2007t|Danilovic et al.|2010HLites|2011| l? These are 
important questions not only in solar physics, but also in stellar physics because there are 
many low-activity dwarf stars and giants with emission line fluxes close to the basal flux 
limit, and in fact define it (e.g. [Perez Martinez et al.|201 l[|Doyle et al.|1994[ l. 

ALMA has the potential to substantially contribute to answering these fundamental 
long-standing questions. By probing the 3D thermal structure and dynamics of the solar 
chromosphere, observations with ALMA give important clues on the sources and sinks of 
heating in the chromosphere and with it essential clues on the heating of the corona since the 
energy required for heating the corona must pass through the chromosphere. Atmospheric 
heating is addressed in more detail in the context of solar flares in Sect. |4.4| and for regions 
with strong magnetic fields in Sect. |4.5[ Chromospheric waves and oscillations, which are 
potentially important sources of heating, are discussed in Sect.|4.6| 
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Solar flares are violent eruptions during which large amounts of energy are released in the 
form of radiation and high-energy particles. Major flares occur in active regions with strong 
magnetic fields. The strongest flares (X class) can - often in connection with coronal mass 
ejections (CME) - directly influence heliosphere and Earth (see Sect. |4.1| ). However, there is 
a large spectrum of flare energies, which reaches down to so-called micro- and nano-flares. 
Such weaker events can occur on smaller spatial scales and potentially are candidates for 
heating the “quiet” solar corona even outside active regions. Despite their obvious impor¬ 
tance, the basic physics of flares at all scales is not understood, especially regarding particle 
acceleration mechanisms, but also concerning the source of the still enigmatic emission 
component at sub-THz frequencies. Furthermore, quasi-periodic pulsations are common 
and possibly intrinsic features of flares. Understanding the nature of the pulsations puts 
important constrains on the interpretation of physical mechanisms responsible for the ac¬ 
cumulation and release of magnetic energy and acceleration of charged particles in flares, 
and may become an important diagnostic tool for the physical conditions in the flaring sites. 
The progress in this research direction requires spatial resolution superior than provided 
by available observational tools, as well as a high spectral resolution in the microwave and 
sub-THz bands. ALMA’s spectral, spatial and temporal resolution for observations in the 
sub-THz range together with the ability to probe the thermal structure and magnetic fields in 
flaring regions promise therefore ground-breaking discoveries in this respect. See Sect. |4.4| 
for more details. 

Solar prominences are regions of cool and dense plasma lying in a much hotter and rarer 
solar corona. While active prominences exhibit large-scale motions and erupt in a matter 
of hours, quiescent prominences have large-scale structures that remain stable for days or 
weeks, although their fine structures change on time-scales of a few minutes. The existence 
of prominences is mainly due to the coronal magnetic fields that support prominence ma¬ 
terial against gravity and insulate it from the hot coronal environment. Prominences are 
referred to as (dark) filaments when observed in absorption against the (bright) solar disk. 

ALMA will provide us with unprecedented details on the temperature structure of the 
prominence material, which has direct implications for the energy balance of prominences. 
ALMA’s high spatial and temporal resolution will also be key in the understanding of the 
small-scale dynamics of the prominence fine structures and to the understanding of the origin 
of prominences. ALMA will provide invaluable insight into the role and impact of promi¬ 
nences on the evolution of the physical conditions in our space environment. See Sect. |4.8| 
for more details. 

“Space Weather" is a very dynamic area of research, which, broadly speaking, deals with 
the study of perturbations of the heliospheric environment, including around Earth, the plan¬ 
ets and other solar-system bodies, originating from the Sun. ALMA observations addressing 
the aforementioned major research questions will have implications for the underlying so- 

prominence (a.k.a. filament) 
eruptions (Sect. |4.8| l, and CMEs, including their triggers and internal structure right from 
the events’ birth. 


lar drivers of space weather, including solar flares (Sect. \AA\ 


4.2 Quiet Sun regions 

Our understanding of the chromosphere in quiet Sun regions has been improved substan¬ 
tially during the last decades owing to advances in observation techniques at ultraviolet 
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(UV), visible and near-infrared (NIR) wavelengths and advances in numerical simulations. 
The cartoon in Fig. [^illustrates the complex structure of quiet Sun regions with a wealth 
of dynamic physical processes as it is known today. Nevertheless, quiet Sun internetwork 
regions are reasonably well understood so that continuum intensity observations of such 
regions may be used for radiometric calibrations of ALMA, which could be of potential 
interest for radio astronomers in general. 

4.2.1 The thermal structure and dynamics of Quiet Sun regions 

The classical semi-empirical reference atmospheres by |Vernazza et al.| ( [1981| hereafter re¬ 
ferred to as VAL) and [Fontenla et al.| l [1993[ hereafter FAL) were constructed by adjusting 
a (hydrostatic) plane-parallel atmospheric stratification until the corresponding emergent 
spectra optimally matched a set of selected observations (see also |Avrett & Loeser|2008| l. 
Among the selected target observations were EUV continuum intensities and the Lyo- and 
Lyf spectral profiles, although the spatial resolution of those observations was rather low 
compared to today’s standards. While elaborate at the time, especially in capturing the NLTE 
subtleties of the optically-thick line formation, this class of static models cannot account for 
the pronounced small-structure and dynamics of the solar atmosphere so impressively docu¬ 
mented by modern high-resolution instruments. The hope of the VAL and FAL-type models 
was that the assumption of a one-dimensional, plane-parallel atmospheric stratification and 
the usage of temporally and spatially averaged observations, would lead to a stratification 
that reproduced the atmospheric properties in the mean. However, the true structuring and 
dynamical properties of the atmosphere, especially in the chromosphere, are so severe that 
an average model cannot capture a meaningful picture of the actual complexity. For exam¬ 
ple, the observation of dark carbon monoxide (CO) lines at the extreme limb implied the 
presence of cold gas, 1000 K cooler than the minimum temperatures of the VAL- and FAL- 
type models (e.g. [Noyes & Hall|1972[[^res & Testerman| 19811 1. In fact, the EUV and UV 
intensities that form the basis of the ID models in the chromosphere have a highly non-linear 
relation to the (local) gas temperature, which renders the average stratifications essentially 
meaningless in the presence of unresolved significant temperature variations. Now, however, 
ALMA provides an unprecedented measure of the small-scale thermal structure of the chro¬ 
mosphere and at the same time potentially exploiting the CO molecules in a new light as 
an additional thermal probe (see also Sect. |2.3^ . Using IRIS, |Schmit et al!] ( |2014| l observed 
absorption features which indicate temperatures below 5000 K existing above the transition 
region. The authors suggested that this lends credence to numerical studies that pointed to¬ 
wards elevated pockets of cool gas in the chromosphere - a hypothesis that can be tested 
with ALMA. 

An early attempt to account for atmospheric dynamics was the ID time-dependent simu¬ 
lations of Carlsson & Stein I jCarlsson & Stein|1994||19^ of chromospheric Call excitation 
by propagating shock waves. These waves were generated from initial disturbances in the 
photosphere, which then propagate upwards into layers with lower density, and ultimately 
steepen into shock waves in the chromosphere. The resulting strong variations in quantities, 
especially the gas temperature, can successfully reproduce spatially and temporally aver¬ 
aged near-UV observations (of the Ca II “flashes” or “K grains”) and at the same time allow 
for cooler, less disturbed intervals in-between shocks, as implied by observations of carbon 
monoxide in the low chromosphere ( [Asensio Ramos et al.|2003| l. The result emphasised the 
ambiguities introduced by relying solely on UV diagnostics to empirically deduce proper¬ 
ties of the solar chromosphere, again owing to the highly non-linear response of the UV to 
local conditions in a thermally diverse outer atmosphere. 
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Fig- 9 Schematic structure of Quiet Sun regions. This illustration is based on the cartoon by [Wedemeyer- 
|B6hm et al.H2009) but has been updated and now contains vortex flows and torsional motions. The horizon¬ 
tally aligned coloured areas represent the layers which are mapped by ALMA. 


The calculations of the emergent thermal free-free radiation at submillimeter and mil¬ 
limeter wavelengths for the |Carlsson & Stem| ( |1994| l simulations produced brightness tem¬ 
peratures that varied strongly in time due to the propagating shock waves and oscillation 
modes, and thus provide a much more reliable proxy for the local gas temperature in the 
continuum forming layers ( [Loukitcheva et al.|2004[ l. The resulting synthetic (sub)millimeter 
maps were compared to interferometric observations of quiet Sun regions by the Berkeley- 
Illinois-Maryland Array (BIMA) at a wavelength of 3.5 mm ( [Loukitcheva et al.|2006||2008) , 
although the pattern seen on the smallest spatial scales in the simulations were not resolv¬ 
able in these observations (see also Fig. for examples of BIMA images). The dynamic 
picture was further supported by studies based on 3D models (e.g.,| Wedemeyer-Bohm et ah] 
|2007[ [Loukitcheva et al.|2015a| l, which in addition exhibited an intermittent and fast evolv¬ 
ing pattern in the synthetic millimeter continuum intensity maps with structures on scales 
down to 0.1 arcsec and time scales on the order of a few 10 s. This small-scale pattern cor¬ 
relates closely with the thermal pattern at chromospheric heights as it was already known 
from earlier simulations i Skartlien et al.|2000{|Wedemeyer et al.|2003b|a||2004[[Hansteen^ 


|Gudiksen|2005HCarlsson et al.|2010' i. The pattern is produced by propagating shock waves 

like those seen in ID simulations, although, crucially, in 3D the waves travel in all directions, 
not just vertically. The dynamic co-existence of hot gas in shock fronts and cool gas in post¬ 
shock regions finally provides a way to understand the presence of cold carbon monoxide in 
the chromosphere, which otherwise would be completely at odds with the high temperatures 
in VAL-type models ( ^res^&RabnJl^% [Asensio Ramos et al.|2003^ j Wedemey er-B dhm| 


et al.|2003{|Wedemeyer-Bbhm & Steffen|2007 1 . 

However, because the molecular formation rates are very slow in the post-shock gas, due 
to the low densities, even considering multi-step chemical pathways ( [Ayres & Rabin|1996^ , 
and because the post-shock regions are expected to be spatially intermittent, it is somewhat 
puzzling that the off-limb emissions of the CO fundamental rovibrational bands at 5 yum are 
so strong and spatially prevalent at the limb. These emissions, originally noted by [SolanM| 
|et al.[ ( |T994| , occur because cold gas apparently exists in the chromosphere. It can appear in 
emission above the limb despite the low temperatures of the gas (~3600 K), because ther¬ 
mally emitting molecules are seen against the very dark IR sky at heights above the limb 
where the 5 jj.m photospheric continuum has become transparent. [Ayres] ([2002[l carried out 
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over-the-limb 2-D spherical transport simulations of the CO 5 /rm rovibrational spectrum in 
a variety of models including: (1) the VAL laterally homogeneous reference model; (2) mul¬ 
tiple slices from the ID time-dependent simulation by |Carlsson & Stem| ( |1994| l, with the ID 
snapshots arrayed at random across the surface with a lateral dimension compatible with the 
dynamic “K-grains” observed prolifically in the quiet Sun (~lMm); and (3) various spa¬ 
tial combinations of the warm VAL thermal profile with cold temperature components in 
the low chromosphere region up to about 1000 km, where the hot chromospheric canopy 
naturally destroys the molecules. In these simulations, the homogenous VAL temperature 
profile produces only very minor limb extensions of the CO bands, coming entirely from 
the Lmin region of that model. The [Carlsson & Stein| ( |1994| l variant did much better, be¬ 
cause of cold regions in the higher altitude range due to adiabatic cooling associated with 
the acoustic pulses, producing significant CO limb emissions, although not as high above 
the limb as seen in observations. Models with VAL-like temperature profiles interspersed 
with small-scale cold regions, with temperatures below 3000 K for most of the 500-1000 km 
height range but with only 20 % covering fraction, did about as well as the |Carlsson & SteTn| 
l |1994[ l model, but still did not reproduce the full spatial extension of the observed strong CO 
lines above the limb. Only a model in which the hot part of the VAL in the chromosphere 
was reduced to 20 % coverage, with the remaining 80 % allocated to a temperature profile 
with a cool dip in the chromospheric layers, although with minimum temperatures above 
3000 K, was able to fully reproduce the details of the CO off-limb emissions observed under 
conditions of very good seeing. These simulations seem to suggest that molecule-friendly 
conditions are pervasive in the low chromosphere, with hot gas being a relatively minor 
constituent. 

The simulations, moreover, predict a very different appearance of the CO lines when 
observed on the disk, say in a heliogram taken in the core of a strong 5 fim CO line. The 
homogeneous VAL model of course would show a smooth featureless CO image across the 
surface, aside from large-scale slowly varying center-to-limb effects. The dynamic model by 
[Carlsson & Stein| ( |l994| l predicts a pattern of light and dark features across the surface, about 
equally mixed, with Doppler shifting of the CO lines in each coherent patch. Near the limb, 
however, the intermittent dark patches, tied together by the nearly horizontal light paths, 
dominate the scene because they are optically thicker in CO than the warmer regions. For the 
two VAL -H COOLC (cool component) simulations, the 20 % COOLC version delivers a disk 
center image that is dominated by the brighter emission from the warm VAL component, 
but has conspicuous dark points from the intermittent but very cold regions, which become 
optically thick in CO in the cold “COOLCl” lower chromosphere. On the other hand, the 
80 % “COOLCO” version shows a more bland image, brighter regions sprinkled between 
larger dark regions. The latter, however, are not as dark as for the minority COOLCl case 
because the less extreme thermal profile of COOLO has the disk-center CO lines becoming 
optically thick near the classical Tmin, where the gas is warmer, rather than higher up where 
the very cool temperatures occur. 

These models admittedly are quite simplistic compared to the 3D dynamical simula¬ 
tions of these complex layers that have become available in recent years (see Sect [4.2.2) . 
Nevertheless, the two-component VAL -i- COOLC simulations of on-disk and off-limb CO 
are instructive, and have important implications for ALMA. Namely, both types of models 
(minority cool mixed with majority hot, and vice versa) predict the off-limb extensions of 
CO, and the behavior close to the limb, with varying degrees of success, but much better 
than traditional ID reference models. However, the two general classes of two-component 
simulations make very different predictions for a disk-center quiet-Sun CO-core heliogram. 
Ironically, perhaps, the minority cool component has the most obvious contrasty appearance. 
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with small-scale very dark spots superimposed on a warmer background; whereas the major¬ 
ity cool component heliogram would appear much blander, with a pervasive, somewhat dark 
background with a sprinkling of small-scale not-so-bright points. In short, the cool regions 
are more visible when they are less important. In principle, contemporary CO imaging could 
settle the question straightforwardly, and indeed the best quality CO heliograms available 
do not reveal any obvious small-scale very dark points needed in the situation where the 
cool component is in the minority. However, that conclusion must be tempered by the fact 
that the spatial resolution of existing solar telescopes at 5 /rm is only somewhat less than 
1 Mm so that smaller-scale very dark points could have easily remained unresolved and thus 
undetected. The 4-m DKIST solar telescope will help resolve the question when it comes 
on line later this decade, but its resolution at 5 fim only will be less than a factor of 3, or so, 
better than existing 1.5-m-class solar instruments. 

On the other hand, ALMA not only can image the complex layers of the low chromo¬ 
sphere in the sub-mm region with a spatial resolution much better than today’s largest solar 
instruments in the thermal-IR, but also tune through those layers by observing at different 
frequencies. A single high spatial resolution time series of surface maps by ALMA at its 
highest frequencies will immediately decide whether cool gas in the lower chromosphere is 
pervasive, or not; and whether its existence is owed to intermittent dynamical phenomena, 
or instead perhaps to more stable pockets of gas that are not strongly heated for a time and 
suffer a “molecular cooling catastrophe”, which plunges their internal temperatures to very 
low values (|Ayres|1981[l. 


Table 1 Properties of the considered numerical models and the resulting synthetic brightness temperature 
maps. Mode l A has been compu ted with the MHD code Bifrost jGudiksen et al.|201 Ij and models B-E with 
CO^BOLD jpreytag et al.|2012} . The models have already been described in the following publications: L15 
Loukitcheva et al.|2015aj , W07 jWedemeyer-Bohm et al.|2007j , W12 jWedemeyer-Bohm et al.|2012j , W13 
Wedemeyer et al.|2013b| . 

The values in parentheses for model A refer to the region outside the circle drawn in Fig.llOji. 


Model_I_A_I B I C I D I E 


Reference 

L15 

W12 

W12, W13 

W07 

W07 

MHD code 

Bifrost 

CO^BOLD 

CO^BOLD 

CO^BOLD 

CO^BOLD 

Radiative transfer code 

L15 

Linfor3D 

Linfor3D 

RH 

RH 

Horizontal cells 

504x504 

286x286 

286x286 

140 X 140 

140 X 140 

Horizontal extent [Mm] 

24.1x24.1 

8.0 X 8.0 

8.0 X 8.0 

5.6 X 5.6 

5.6 X 5.6 

Horizontal extent [arcsec] 
Initial / avg. photospheric 

33.3x33.3 

ll.Ox 11.0 

ll.Ox 11.0 

7.7 X 7.7 

7.7 X 7.7 

magnetic field strength 
Non-equilibrium H 

50 

50 

0 

0 

0 

ionisation 

yes 

no 

no 

yes 

(yes) 


Synthetic brightness temperature at i = 1 mm 


Average [K] 

4785.9 (4673.6) 

4568.6 

5086.9 

3935.9 

4476.5 

Standard deviation [K] 

822.9 (736.2) 

879.7 

1085.5 

934.6 

998.2 

Minimum [K] 

2705.8 (3158.4) 

2966.6 

3044.2 

2686.7 

3109.9 

Maximum [K] 

11093.8(8928.7) 

7481.3 

7257.4 

7272.8 

7081.7 

Contrast 

17.2% (15.8%) 

19.3 % 

21.3% 

23.7 % 

22.3 % 
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Fig. 10 Synthetic brightness temperature at i = 1 mm for model snapshots from four different simula¬ 
tions (see Table From top to bottom: A Bifrost simulation with an average magnetic field strength of 
< |B| >ph = 50G in the photosphere (model A), and CO^BOLD simulations with an initial field strength of 
|B|o = 50G (model B), the corresponding non-magnetic model C, and the smaller model D, which takes into 
account time-dependent hydrogen ionisation. The upper left panel shows the primary beam size of an 12-m 
ALMA antenna (i.e. the diffraction limit), which has a diameter of 21” at /I = 1 mm (dot-dashed circle) and 
the sizes of the smaller CO^BOLD models (dashed rectangle for models B and C and dotted rectangle for 
model D). The corresponding brightness temperature histograms are shown in the right column as light grey 
area. For comparison, a second distribution is displayed as a dark grey area. For model A (panel e), the his¬ 
togram of the pixels inside the primary beam (the circle in panel a, dark grey here) are compared to those 
outside (light grey), where the latter excludes the bright loops in the center. The distribution for the whole 
model is over-plotted as red thick line. In the panels f and g, the models B and C are compared to each other. 
In panel h, model D (light grey, NLTE) is compared to model E (dark grey, NLTE), which is equivalent to 
model D but uses LTE electron densities instead of the results from the time-dependent hydrogen ionisation. 
The vertical lines mark the mean brightness temperature, whereas the horizontal lines with square symbols 
at the top of each panel represent the mean plus/minus one standard deviation. The histogram for the pixels 
outside the circle in model A (i.e. without the magnetic loops in the center) is plotted as reference in panels f-h 
(red triple-dot-dashed line). 
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4.2.2 Numerical predictions of quiet Sun ALMA observations 


In Figure [To] and Table [T] we compare different 3D numerical models. Model A was com¬ 
puted with the Bifrost code ( [Gudiksen et al.|201 l[ l and is classifie d as ‘enhanced netw ork’ 
whereas the other models (B-D) were calculated with CO^BOLD ( [Frey tag et al.|2Q12 i and 
represent quiet Sun regions. The synthetic brightness was produced with the rad iative trans- 
fer codes LINFOR3D (based on the Kiel code LINFOR/LINLTI^ and RFl ( Uitenbroek 
|2000 1 , whereas the results for model A were calculated in detail by |Loukitcheva et al7 
1 2015a| l. Model A has a horizontal extent of 24.1 Mm x 24.1 Mm (33”.3 x 33”.3) and reaches 
into the corona. The magnetic field strength in the photosphere is on average < |B| >ph = 
50 G with two major regions of opposite polarity, which are roughly 8 Mm apart. The photo- 
spheric patches of opposite polarity are connected by magnetic loops, which extend into the 
low corona. These magnetic loops are nevertheless visible in the corresponding synthetic 
brightness temperature map for a wavelength of T = 1 mm (see upper left panel in Fig.|10|i. 
The brightness temperature in these loops is higher than in the surrounding (quiet) chro¬ 
mosphere but only 0.3 % of all pixels have a brightness temperature higher than 8000 K. 
The corresponding histograms are displayed in the lower left panel for the whole simula¬ 
tion and the outer four corners, which thus exclude most of the bright loops in the centre. 
The brightness temperatures range from only ~ 2700 K to ~ lllOOK with an average of 
~ 4790 K and a standard deviation of ~ 820 K. The 1 mm map is compared to a corre¬ 
sponding 3 mm map and Mg Ilk line wing and line core images in Fig. The numerical 
model A is described in more detail in |Carlsson et al.| l f2016[ see also IRIS technical note 33, 
http://sdc.uio.no/search/simulations). 

The other models have a smaller horizontal and vertical extent and therefore do not 
include large-scale magnetic loops like those in model A. The maximum brightness temper¬ 
ature at /I = 1 mm is therefore 7 500 K or less in models B-D without the high-temperature 
tail in the histograms, which is caused by the hotter loops in model A. Models B, C, and D 
nevertheless show differences due to different magnetic field, spatial extent and inclusion 
of non-equilibrium hydrogen ionisation. Models B l [Wedemeyer-Bdhm et al.|2012t |Wede-| 
|meyer et al.|2013b| and C l [Freytag et al.|201^ have the same spatial extent but model B 
has an initial magnetic field strength of |B|o = 50G whereas model C is non-magnetic. The 
histograms of the brightness temperature at T = 1 mm are over-plotted for both models in 
the middle bottom panels of Fig.[T^so that the influence of magnetic fields becomes direcfly 
visible. The pronounced small-scale mesh-like patfern, which is produced by the interaction 
of propagating shock waves in model C (see third brightness temperature map in panel c), 
is somewhat suppressed by the magnetic field structures in model B. These field structures 
are rooted in kilo-Gauss foot points in the photosphere and expand into chromospheric fun¬ 
nels, which then influence the propagation of waves in the upper layers. Consequently, the 
peak in the histograms at high brightness temperatures due to shock fronts is pronounced for 
model C but ablated in model B. On the other hand, chromospheric swirls appear in the map 
for model B (panel b) and indicate the presence of magnetic tornadoes (see Sect. |4.2.4| l. 

In Fig. 1 1 op, the corresponding brightness temperature for the non-magnetic model D is 
shown, which has the smallest horizontal extent of the models considered here. For model D, 
the hydrogen ionisation is treated in detail by solving the time-dependent rate equations for 
a six level model atom with fixed radiative rates l [Leenaarts & Wedemeyer-Bohm|2006| l. The 
resulting non-equilibrium (non-LTE) electron density varies much less in the chromosphere 
than the corresponding LTE electron densities, which instantaneously follow any change 


LINFOR3D manual at http://www.aip.de/~mst/linfor3D_main.html 
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Fig. 11 Spatial resolution effect on the intensity histogram for the wavelength of 1 mm and four values of 
FWHM (spatial resolution): 0.”2 (solid), 0.”4 (dotted), 1” (dashed), and 4” (dot-dashed). 


in local gas temperature. The finite recombination rates in the non-equilibrium case effec¬ 
tively weaken the variations caused by the hot propagating shock waves. The brightness 
temperature distribution of the non-LTE case is compared to the corresponding LTE case, 
i.e. the same model but with equilibrium electron densities, which is referred to as model E 
in Fig.[T0]t and Table[T] An analysis of the contribution functions of the continuum intensity 
reveals that the variation of the effective formation height for different spatial locations at a 
given wavelength in model E is reduced to 60 % compared to model D, demonstrating the 
importance of a realistic treatment of the hydrogen ionisation. The contribution functions 
for model E are shown in Fig.[^ According to |Wedemeyer-B6hm et al.| ( |2007[ l, the narrower 
continuum forming layer in model E results in the reduction of the average brightness tem¬ 
perature in the resulting maps from 4,476 K in model D (i.e. with LTE electron densities) 
to 3,935 K for model E (i.e. with the more realistic non-equilibrium electron densities). The 
standard deviation also decreases somewhat from dTb LTE = 998 K to dTb.NLTE = 933 K. 
In both cases, the histograms exhibit a very pronounced peak at low temperatures. The 
high-temperature peak due to shock waves is reduced when non-LTE electron densities are 
considered, suggesting that ALMA will detect a considerable amount of cool gas that goes 
unnoticed by atomic lines in the visible and UV. 

Spatial resolution emerges as a critical problem for observations of the chromospheric 
fine-structure with ALMA. |Loukitcheva et ah] ( |2015a[ l studied the effect of limited spatial 
resolution by applying artificial smearing to synthetic millimeter brightness temperature 
maps from a 3D snapshot of the aforementioned model A. The minimum resolution in the 
model is defined by the horizontal extent of the grid cells (0.”064, cf. Table[T|. |Loukitcheva| 
|et al.[ found that the fine structure on the smallest spatial scales in the synthetic millimeter 
maps is increasingly lost when the resolution is reduced. At a resolution of 1” the brightness 
temperature contrast decreased to 70-80 % of the original contrast value, depending on the 
wavelength. The fine structure on the smallest scales is already lost, while the larger-scale 
pattern is still distinguishable and closely reflects that in the original image. In the intensity 
histograms shown in Fig. [m these effects are visible as narrowing of the histograms, disap¬ 
pearance of the high-brightness tail and growth of the low-brightness bulk. These findings 
are in line with the earlier study by|Wedemeyer-Bohm et al.|(|2007| see their Fig. 10). They 
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Fig. 12 Longitudinal component of the magnetic field at the photospheric level (left) and as it is derived 
from synthetic millimeter brightness temperature maps at the chromospheric level (right). The latter roughly 
corresponds to a height of 1500 km, which is where the emission at d = 3 mm mostly comes from in this case. 
The value range for photospheric magnetic field is limited to [-500 G, -t500 G] for clarity. The employed 
brightness temperature maps are calculated based on a snapshot of the 3D RMHD simulations done with the 
Bifrost code (Gudiksen et al.|201 1) . 


State that the pattern on the smallest scales are hardly discernible at a resolution lower than 
0.”3 whereas the larger mesh-like pattern due to shock waves can still be seen at a resolu¬ 
tion of 0.”9. |Wedemeyer-Bdhm et SI] approximate the reduction of brightness temperature 
contrast 6T^,^nns/{Th} as function of angular resolution Aa with an exponential of the form 

^ (20) 

(^b) 

The characteristic length scale D of the chromospheric pattern is derived from fits to the 
contrast curves and depends on wavelength and view angle. For a wavelength of 1 mm at 
disk centre, for instance, the contrast reduces from 24 % to just 10 %, i.e. to about 40 % of 
the original value, when imposing artificial smearing with a spatial resolution corresponding 
to 1”. 


4.2.3 Magnetic fields in Quiet Sun regions 


|Loukitcheva et al.| ( |2015b| ) carried out tests of the technique described in Sect. 2.5 for ALMA 
frequencies by using a 3D simulation snapshot from the Bifrost code (Gudiksen et al.|201 1| ) 
like the one described in Sect. |4.2.2[ For a region of enhanced network with maximum 
magnetic field of 2kG at the photosphere (left panel in Fig.|12|l and its rapid decrease with 
the height, the simulated circular polarisation at 3 mm falls within ±0.5 %, which results in 
the restored longitudinal magnetic field within ± 100 G at the effective formation height. The 
(expected) polarisation is thus higher than the technical requirement of 0.1 % for ALMA 
(see Sect. |3.3.4[ l and could thus be used to measure the magnetic field in the quiet Sun 
chromosphere. 

An example of the magnetic field strength at the height of approximately 1500 km, as it 
is restored from the simulated polarised emission at T = 3 mm, is shown in the right panel of 
Fig. For this example, the difference between the effective formation heights of the two 
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normal modes does not exceed 2 km at d = 3 mm and is only slightly greater than 3 km for 
A = 10 mm. The assumption of a constant magnetic held, which is implicitly made for this 
method, is valid for such small separations. Consequently, the method reproduces the longi¬ 
tudinal magnetic held with the ideal model data reasonably well. However, further analysis, 
e.g. involving realistic noise, is needed to fully validate the approach. Corresponding results 
for predicting magnetic helds in active regions are discussed in Sect. |4.3.T] 

4.2.4 Vortex flows 


A magnetic vortex how is a fundamental physical process in the solar atmosphere, which 
forms when the photospheric footpoint of a magnetic held structure is swept into a vortex 
flow in an intergranular lanes at the solar surface (see Fig. 0. The primary photospheric 
vortex is a direct consequence of the conservation of angular momentum of the downhow- 
ing plasma in the intergranular lanes (‘bathtub effect’) and thus an integral part of stellar 
surface convection (see, e.g., |Bonet et al.||2008[ |Kitiashvili et al.||2011| |Moll et al.||2011| 
|2012[ and references therein). The rotation of the magnetic footpoint is mediated into the 
atmosphere above, where the plasma is forced to follow the rotating held structure (|Wede-| 
[meyer & Steiner|| 2014[l. There, it forms a observable sig nature, namely a ‘chromospheric 
swirl’ jWedemeyer-Bohm & Rouppe van der Voort|2009t This kind of event, which is also 
referred to as a ‘magnetic tornado’, may serve as a channel for energy and matter into the 
chromosphere and corona of the Sun and has therefore a potential contribution to heating 
the upper solar atmosphere ( |Wedemeyer-Bohm et al.|2012[ l. Based on measurements with 
the Swedish 1-m Solar Telescope, chromospheric swirls have lifetimes of (12.7 ± 4.0) min 
and diameters of (4.0 ± 1.4 arcsec (i.e, (2900 ± 1000) km), whereas the width of substructure 
features (e.g., rings or spiral arms) can be as small as < 0.”2. ALMA observations of swirls 
(cf. Fig.|10[?) would provide a direct measure of the gas temperatures within these structures 
and thus a mean to further investigate their role for the chromospheric energy transport and 
the possibly resulting heating of the outer solar layers. 

Rotating magnetic field strucfures are observed on an extended range of different spatial 
scales. On the Sun, it ranges from torsional motions in spicules (|De Pontieu et al.|2012 


to 


very large vortices, which appear in connection with solar prominences (e.g. Li et al.|2012^ . 
It is currently debated if the “legs” of prominences, also referred to as ‘giant solar tornadoes’, 
are rotating, too (e.g., |Wedemeyer et al.|2013a| l. See Sect. |4.8| for more information on solar 
prominences and Sect. 4.9 for potential implications for other helds of astrophysics. 


4.2.5 Modelling of ion-neutral collisions 


State-of-the-art numerical simulations so far employ a single-fluid (ideal) MHD description, 
which is strictly speaking only valid if the relevant time scales are larger than the ion-neutral 
collision time. The latter is very short in the photosphere (10“® - 10“® s) and the lower chro¬ 
mosphere (10“^ - 10“^ s). Therefore, in this part of the solar atmosphere the single-fluid 
description perfectly describes the processes in the limits of spatial and temporal resolutions 
of ALMA. In middle and higher chromosphere, however, the chromospheric plasma is par¬ 
tially ionised with the ionisation degree varying significantly l |Carlsson & Stein|2002[ l. The 
collision time between ions and neutral atoms (especially, neutral helium) is dramatically in¬ 
creased in the upper chromosphere and may reach the scales of ALMA observations (~ 1 s 
for proton-neutral helium collision). Therefore, the validity of the single-fluid approach has 
to be critically checked through comparisons of simulation results with ALMA observations. 
Such studies may motivate moving forward towards a multi-fluid MHD approach (see, e.g. 
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|Zaqarashvili et al.|2011[ l, which is ultimately wanted for a realistic modelling of the propa¬ 
gation and dissipation of Alfven waves in the solar atmosphere (see Sect. |4.6.^ . 

4.2.6 Center-to-limb variation and spicules 


Observations in the mm and cm range show less limb brightening than expected from 
the increase of the temperature with height ( [Lantos & Kundu||1972[ [Bastian et al.| | 199?b| 
|Kochanov et al.|2013| l and this cannot be interpreted in terms of the instrumental resolution. 
It is usually explained as the effect of absorbing features, and spicules have been implicated 
to this end. In order for spicules to absorb, they have to be cooler than the ambient medium, 
which is certainly true for the transition region and the corona into which they penetrate, but 
not necessarily so in the chromosphere (see, e.g., Tsiropoula et al.|2012[ l. 

In the 2D atmospheric model proposed by |Selhorst et al. (2005a i, the authors included 
the presence of spicules to explain the low equatorial limb brightening (~ 10%) observed at 
17 GHz. On the other hand, the high values observed at the poles were modeled by spicules- 
free regions associated with polar faculae ( |Selhorst et al.|2005b[ l. 

In order to measure the center-to-limb variation with ALMA the TP system must be fully 
implemented; then one can measure the brightness of both the network structures and the 
background at various positions on the disk and at the limb. Even without TP, the brightness 
difference between network and intranetwork can be measured and this is an important di¬ 
agnostic. Spicules themselves and the corresponding disk structures will be hard to observe, 
because their short lifetime (~ 15 min for normal spicules, ~ 2 min for type II spicules; see 
[Tsiropoula et al.|2012[ l will probably not allow sufficiently dense coverage of the u-v plane 
to properly image the complex structure with the bright limb in the field of view. 


4.2.7 Polar brightenings 


Radio observations have revealed the remarkable result that the polar regions of the Sun are 
brighter than the rest of the quiet Sun in a limited range of frequencies from 17 GHz to about 


87 GHz (e.g., Efanov et al.|1980[|Kosugi et al 

|1986 

|Shibasaki|1998[ Nindos et aL|1999b| 

Pohjolainen et al.||2000 Selhorst et al.||2003| 

2010 

2011 Gopalswamy et al.||2012[ Nitta| 

et al.|2014i. Most of the recent studies of the polar brightening have used data obtained by 


the Nobeyama radioheliograph (NoRH) at 17 GHz. The NoRH images show that the polar 
brightening consists of two components: a diffuse component of 1500 K excess brightness 
and patchy compact sources with localized excess brightness of about 3500 K. The compact 
sources change significantly with time. While some weak short-lived features (lifetimes 
about 1 minute) are difficult to distinguish from side-lobes, some others live longer, do not 
move, and are probably real sources l| Nindos et al.|1999b| l. 

|ShibasakT| ( |1998[ l proposed that the polar brightening arises from a general limb bright¬ 
ening (due to an enhanced opacity in the line of sight) superposed on the bright features 
intrinsic to the polar region. However, the nature of the bright compact sources is largely 
unknown and they have no correspondence with small-scale bright regions in images from 
both the Extreme Ultraviolet Imaging Telescope (EIT) aboard the Solar and Heliospheric 
Observatory (SoHO) ( [Nindos et al.|1999b ]l and the Atmospheric Imaging Assembly (AIA) 
aboard the Solar Dynamics Observatory (SDO) ( |Nitta et al.|2014[ l. 

This lack of one-to-one correspondence between 17 GHz and EUV bright regions could 
be due to the presence of radio optically thick features, like spicules, that are completely 
transparent at EUV lines (|Selhorst et al.|2010} and this shadowing may extend onto the disk 
at the extreme limb (e.g., [Lindsey & Hudson] 1976|l. 
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The long-term variation of the polar brightening has also been studied (| Shibasaki| 1998[ 
[Gopalswamy et al.|2012[|Nitta et al.|2014[ l. At 17 GHz the polar brightening anti-correlates 
with the solar cycle as defined by the sunspot number and correlates with the polar magnetic 
field strength. In contrast, no polar brightening is apparent in the synoptic maps obtained 
by SSRT at 5.7 GHz (Altyntsev 2015, private comm.). [Selhorst et al.| ( |201ll l,|Gopalswamy| 
|et al.| ( (20T2| l and |Nitta et al.| ( (20T4| l have found that both the polar brightening and polar field 
strength were larger during the cycle 22/23 minimum than during the cycle 23/24 minimum. 
Concerning the occurrence of minimum polar brightening, they found that the north and 
south solar hemispheres were not tightly synchronized. 

The study of the solar cycle dependence of the polar brightening is only meaningful us¬ 
ing data from solar-dedicated instruments. However, ALMA with its unprecedented spatial 
resolution and sensitivity has the potential to provide important new observations that will 
help us understanding the nature of polar brightening’s compact sources. 


4.3 Active regions and sunspots 

4.3.1 Active region modelling and predictions for ALMA 


As has been detailed in Sect. |4.2[ the chromosphere is extremely difficult to model in the 
quiet Sun, let alone in active regions where the magnetic field becomes even stronger and 
more highly structured. Detailed modeling of the microwave emission from active regions 
(see, e.g.,|Alissandrakis et al.|19^0||Wang et al.|2015 and references therein) showed the 


importance of the gyroresonance mechanism (see Sect. 2.4^ and its potential role as an accu¬ 
rate diagnostic of the magnetic field. In the ALMA spectral range, however, the gyroresonant 
emission is unimportant in the vast majority of the cases, while the free-free emission makes 
the dominant contribution (see Sect.j^for details). 

To support solar ALMA observations of active regions, [Fleishman et al.| ( |2015|l devel- 
oped an efficient algorithm integrated in the 3D modeling tool GX Simulator ( Nita et ah] 
|2015[ available in SolarSoft), facilitating the quick computation of synthetic maps in the 
ALMA spectral range which give a rough idea of how a given active region will look at 
(sub-)millimeter wavelengths. Based on photospheric white-light images and magnetograms 
of a selected active region, an appropriate ID model atmosphere is chosen for each pixel, 
i.e., for each spatial position in the image. The algorithm employs the set of ID static model 
atmospheres by |Fontenla et al.| ( |2009[ |2014[ l, which describe averaged atmospheric strati¬ 
fications corresponding to different types of features in the photosphere including umbra, 
penumbra, network, internetwork, enhanced network, facula, or plage. Each model provides 
the height stratification of all relevant physical quantities including ID distributions of the 
electron temperature and density together with (non-LTE) ionized and neutral hydrogen den¬ 
sities. After the appropriate model is chosen, the emergent millimeter radiation is computed 
for each pixel accordingly. The algorithm accurately takes into account gyroresonance (GR) 
and free-free processes, including magnetized plasma dispersion and frequency-dependent 
mode coupling ( [Fleishman & Kuznetsov|2014| l, non-LTE hydrogen ionization, all heavier 
singly-ionized atoms, the helium ionization state, and emission due to collisions of electrons 
with neutral hydrogen and helium (H“ and He“ opacities). 

Some results of this modeling for AR 12158 (observed 10 Sep 2014) are shown in 
Fig.[T^for two frequencies representing the currently operating ALMA bands 3 and 6. This 
modeling suggests that different features will be distinguishable at the ALMA frequencies 
with the brightness temperature contrast up to a few hundred K. Interestingly, in line with 
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Fig. 13 Examples of simulated images for AR 12158 in two ALMA frequency bands (identified at the top-left 
of the leftmost images). From left to right the plots display brightness temperature in right circularly polar¬ 
ized light, brightness temperature in left circularly polarized light and the relative difference between the two 
polarisations. Note that the umbra appears bright (hot) at lower frequencies but dark (cool) at higher frequen¬ 
cies. The degree of polarisation along with the spatially resolved brightness temperature spectra is the key to 
obtaining the magnetic field map at a given frequency, which translates to a certain (range of) chromospheric 
heights. Multi-frequency data offer an elegant chromospheric tomography method (see, e.g., |Wedemeyer-| 
|B6hm et al.|2007||Loukitcheva et al.|2014||2Q15b) . The QR codes provide links to on-line movies. 


the results of |Loukitcheva et ar] l |2014| t, the umbra, which is dark (cool) at high frequencies 
as well as in white light, appears bright (hot) at lower frequencies. In addition, the bright¬ 
ness temperature noticeably depends on the magnetic field value due to the dependence of 
the free-free opacity on the magnetic field, which can be exploited for magnetic field mea¬ 
surements (see Sect. [g. Obviously, this simplified modeling is only a first step towards 
more realistic, eventually 3D models. In particular, dynamic small-scale processes, which 
can have significant impact on the atmospheric structure, can by nature not be described 
with static ID models (see also the similar discussion for quiet Sun regions in Sect. |4.2.T] i. 
On the other hand, a consistent 3D numerical model of a whole active region with sufficient 
spatial resolution is computationally not feasible with the resources available today. An al¬ 
ternative approximative approach would be to use a systematic grid of smaller 3D MHD 
models, which represent different sub-regions in an AR, ranging from weak-field regions to 
strong magnetic field concentrations. Based on magnetograms of real ARs, like those ob¬ 
tained with SDO/HMI, these different models can then be used as building blocks to fill a 
larger computational domain. 
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In a non-flaring active region the degree of circular polarisation reaches a few percent at 
millimeter wavelengths ( [Fleishman et al.|2015^ , which allows reliable determination of the 
longitudinal component of chromospheric magnetic held from ALMA observations. Exam¬ 
ples of AR polarized emission maps in ALMA frequency bands (left circular polarisation 
(LCP) and right circular polarisation (RCP)) together with the degree of circular polarisation 
(P), reaching 3.6 % at 3 mm, is shown in Fig.[^ For more advanced AR magnetometry a 
combination of ALMA observational data with gyroresonant emission measurements in the 
microwave range available from the Very Large Array (VLA) and, in near future, from the 
Expanded Owens Valley Solar Array (EOVSA), the Siberian Solar Radio Telescope (SSRT), 
and the Chinese Solar Radio Telescope (CSRT), is desirable. The advantage of such a com¬ 
bination is that the gyroresonant emission yields the absolute value of the magnetic field at 
the transition region level (e.g., |Gary & Hurford| 1994^ [Tun et al.|2011[|Wang et al.|2015^ , 
which complements the line-of-sight magnetic field component measurements made with 
the free-free emission. However, co-observing with some of these facilities might be limited 
given the geographic location of ALMA and resulting time differences and, in some cases, 
might be limited to studies of slowly varying features. 


4.3.2 Structure and dynamics of sunspot umbrae 


Although sunspot umbrae are well-studied, there are still a number of important issues that 
have defied attempts at understanding them (see, e.g., the review by |Solanki|2003| l. Among 
these is the rather basic question of the temperature structure of the chromosphere above 
sunspot umbrae and the heating mechanism producing it. Models differ significantly from 
each other, and also the amount of heating required ( jLoukitcheva et al.|2014| l. Studies based 
on data gathered at visible wavelengths give controversial results regarding the height at 
which the chromospheric temperature rise starts in umbrae. [Loukitcheva et al.|p014| l have 
shown that it is possible using millimeter-wavelength observations to distinguish between 
the different models. This approach was so far hampered by the low spatial resolution of 
available mm observations. This problem will be easily overcome by ALMA, even if not 
used in its highest resolution configuration. Furthermore, the ability of ALMA to observe 
at different wavelengths and thus to sample different atmospheric layers will facilitate exact 
measurements of the thermal structure above sunspot umbrae, including determinations of 
the exact height of the temperature minimum layer. 


Sunspot umbrae have turned out to be far more complex than previously thought with a 
wealth of (often hidden) fine structure and dynamical phenomena (see, e.g.,|Rouppe van d^ 


Voort & de la Cruz Rodrfguez|2013[ Yurchyshyn et al.|2015|l, such as umbral dots (|Beck~ 


& faflant|1969t Wittmann|1969t|Riethmuller et al.|2013 i, umbral flashes ( [Beckers & Tallant 


1969[ de la Cruz Rodrigu ez et al. |2013b i, umbral oscillations ( [Sych & Nakariakov|2014^ , 
umbral microjets ([Bharti et al.|2013|, light bridges of different types and their oscillations 


(e.g., Rimmele 1997[ Lagg et al.|2014[[Yuan et al.|2014[l, jets emanating from light bridges 


( ShimiziJ20n i as well as umbral brightenings which are footpoints of hot coronal loops 
( Alissandrakis & Patsourakosj2013[ l. Many of these phenomena (flashes, oscillations, mi¬ 
crojets and light-bridge jets) are clearly chromospheric phenomena whose study will greatly 
profit from the ability of ALMA to trace different temperature regimes and heights as well 
as from the extremely high cadence that ALMA provides. The high cadence allows the evo¬ 
lution of these often short-lived phenomena to be followed in detail. 
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4.3.3 Umbral oscillations and penumbral waves 


Umbral oscillations (which are actually waves, see |Alissandrakis et aL|1992^ were detected 
in the microwave range by |Gelfreikh et al.| ( |1999[ l in Nobeyama circular polarisation data at 
17 GHz; they were spatially resolved in VLA observations both in Stokes I and V at 8.5 GHz 
and 5 GHz ( Nindos et al.|2002) l where they appeared as localized, intermittent oscillations. 
[Reznikova et al.| (2012 1 analyzed a long time series of simultaneous NoRH and AIA obser¬ 
vations and determined time lags and spectral characteristics of umbral oscillations, which 
were consistent with those of upward propagating magnetoacoustic gravity waves, with the 
magnetic held acting as a waveguide up to the ~ 1 MK level. 

Penumbral waves, which may l |Tsiropoula et al.|2000^ or may not ( [Christopoulou et ah] 
|2001| l be an extension of umbral waves have not been detected in the radio range, probably 
due to the very intense sunspot-associated emission and/or the scarcity of high resolution 
observations. Recent observations (e.g. [Madsen et al.|2015l l are consistent with the picture 
that umbral and penumbral waves have a common origin and propagate upwards along held 
lines the inclination of which varies between the sunspot center to the penumbra. 

Both umbral oscillations and penumbral waves are expected to be associated with den¬ 
sity variations that produce changes in the free-free emission detectable by ALMA; having a 
relatively simple spatial structure, they should be readily observable in snapshot mode. Such 
observations will ohFer accurate diagnostics of the physical condition in a range of heights 
that will be extremely useful in modeling their behavior. 


4.3.4 Small-scale dynamic events in sunspot penumbrae 


Modern ground- and space-based telescopes have revealed the presence of many dynamic 
events in and above sunspot penumbrae, namely, running penumbral waves, localized up- 


and downhows, inward motion of bright penumbral grains, chromospheric microjets i Kat-| 
[sukawa et al.||2007^, and moving bri ght dots in the transition region ( |Tian et al.||2014by 
and in the corona ( Alpert et al.|2015y Here we focus mainly on the two latest discoveries 
of small-scale dynamic events in penumbrae: penumbral microjets (PMJs), and penumbral 
bright dots (BDs). Examples of both types of events can be seen in Fig.|14|t-d. 

It was thought earlier that penumbral microjets could be generated by the reconnec¬ 
tion between the two magnetic components of penumbra, namely spines (with more vertical 
field) and interspines (with more horizontal field). Recently, [Tiwari et al.|p015) l proposed 
a modified formation mechanism of penumbral microjets based on the investigation of the 
small-scale structure of penumbral filaments by Tiwari et al. ( 2013| l, who performed a spa¬ 
tially coupled inversion ([van Noort|2012[|van Noort et al.|2013 i of spectropolarimetric data 


of a sunspot observed with the Spectropolarimeter (SP) of the Solar Optical Telescope (SOT, 
[Ichimoto et al.|2008)|Tsuneta et al.|2008[ l onboard the Hinode spacecraft. [Tiwari et al.| ( [2015l ) 
also searched for coronal and transition region signatures of penumbral microjets but hardly 
found such signatures except for some larger penumbral jets, which appeared repeatedly at 
the same locations in chromospheric CallH line images. This difference and the origin of 
both, the ‘normal-sized’ and the larger jets, remains elusive. Based on the results of [Tiwari] 
[et al.[ ( [2015[ l and in line with the earlier proposal by [Magara[ ( |2010| l, it seems reasonable 
to assume that both types of jet are produced by magnetic reconnection events occurring 
between spines and opposite polarity fields, which are observed near the heads of penum¬ 
bral filaments. However, only one third of the filaments studied by [Tiwari et al.[ ( [2013| ) 
displayed opposite polarity fields near their heads. On the other hand, much of the internal 
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Fig. 14 Hi-C observation at a wavelength of 19.3 nm taken on July II, 2012. (a) Full field-of-view. (b) Close- 
up region (blue box in panel a) showing an enlarged part of a sunspot penumbra, (c) Corresponding running 
difference (RD) image of the close-up region showing many penumbral bright dots jAlpert et al.|2015t , which 
also can be identified in panel b. (d) A running difference image in the CallH line obtained by Hinode’s 
filtergraph for the same field of view, exhibiting a few penumbral microjets (seen as white streaks). These jets 
are studied in detail by|Tiwari et al.|{2015). 


Structure of penumbral filaments might have remained undetected due to the limited spa¬ 
tial resolution of the employed instruments. To explore the small-scale internal structure of 
sunspot penumbrae, observations with very high spatial resolution are therefore crucial. In 
this respect, tracking dynamic features at high temporal cadence at (subarcsecond) angu¬ 
lar resolution with ALMA could provide important clues on the formation mechanism and 
evolution of penumbral microjets. Ideally, such observations should be carried out simulta¬ 
neously with observations with other space-borne telescopes such as, e.g., Hinode and the 
proposed Solar-C. 


Penumbral bright dots, which are another kind of dynamic event in sunspot penumbrae, 
move inward or outward along penumbral filaments. Such events were first observed with 
the Interface Region Imaging Spectrograph (IRIS) and then also with the High-resolution 
Coronal Imager (Hi-C, [Kobayashi et al.||20I4| l in a narrow passband range centered at a 
wavelength of 19.3 nm l Alpert et al.||2015 i. See the Hi-C image in Fig. [1^ for some ex¬ 
amples. The fact that the bright dots are detected by IRIS and in the 19.3 nm passband of 
Hi-C, which includes some transition region emission, implies that these events originate 
in the transition region or possibly even lower in the atmosphere, i.e. in the chromosphere. 
The light curves of some of the bright dots observed with Hi-C indeed clearly indicate their 
origin in the transition region, while the light curves for many other events remain incon¬ 
clusive. Among the various possible explanations of the origin of penumbral bright dots 
are, e.g., magnetic reconnection in inclined penumbral magnetic fields, shocks produced by 
downward flowing gas, and reconnection events caused by penumbral waves l [Tian et 3T] 
|2014b[[Alpert et al.|2015|l. Higher cadence and spatial resolution observations with ALMA 
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Fig. 15 SUMER O VI 103.7 nm slit data showing both blue and red shifts events (commonly termed explo¬ 
sive events) over a two minute interval. From|Perez et al.H1999). 


can certainly provide important constraints on the so far elusive origin of penumbral bright 
dots. 


4.3.5 Ellerman Bombs 


Ellerman Bombs (EBs) are usually identified as prominent small-scale brightenings in the 
wings of the Ha line ( |Ellerma^|1917[ l. They have mean lifetimes of ~ 10 - 15 min and 
are mainly observed in the vicinity of active regions ( |Georgoulis et al.|2002^ or regions of 
emerging magnetic flux ( |Watanabe et al.|2008 1 and do not produce any observable effects 
in the transition region or corona ( Vissers et al.|2013| [Nelson et al.|20T5) l. These events are 
generally considered to be due to reconnection in the upper photosphere, which leads to 
increased temperatures therefore explaining the enhanced Ha wings (see e.g. |Rutten et al.| 
|2013| l. High cadence coupled with high spatial resolution data from ALMA, coupled with 
its ability to observe different chromospheric layers can help to answer the question re¬ 
garding the contribution to the chromosphere. For example, do these small-scale, explosive 
phenomena in the lower solar atmosphere act as agents for triggering local instabilities in 
the magnetic environment of the solar surface? 


4.3.6 Explosive Events 


Another feature that ALMA is ideally suited for are the so-called “explosive events” (EEs, 
see Fig.|15[l. These are c haracterized by non-Gaussia n line profiles showing Doppler veloc¬ 
ities of 50 - 150kms“* (Brueckner & Bartoe|l983 i. However, despite many attempts their 
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Fig. 16 Left and middle panel: broadband (1-405 GHz) spectra of a solar flare on 2006 Dec 6 for two dif¬ 
ferent time frames jKaufmann et al.|2009b^ . The data were obtained with the Owens Valley Solar Array 
(OVSA) and the Solar Suhmillimeter Telescope (SST). A mysterious distinct suh-THz spectral component is 
apparent at 212 GHz and 405 GHz. Rightmost panel: broadband (1 -405 GHz) spectra of another solar flare on 
2001 Aug 25 for a few consecutive time frames jRaulin et al.|2004} . In a great contrast with the left panels, 
the radio spectrum of this flare represents a single, presumably, gyrosynchrotron spectrum with unusually 
high spectral peak frequency, around 100 GHz at the peak time, while the spectrum decreases at 212 GHz and 
405 GHz. 


true nature remains unknown, although there is strong evidence that they are produced by 
magnetic reconnection ( |Innes et al.|1997[ l. Recently, |Huang et al.| ( |2014b| using the Interface 
Region Imaging Spectrograph (IRIS) found that EEs are associated with very small-scale 
plasma ejection followed by retraction in the chromosphere. They observed a jet in a single 
IRIS pixel and conclude that it is thus only < 120 km wide. These small-scale jets origi¬ 
nate from a compact bright-point-like structure of ~ l."5 size as seen in the IRIS 133.Onm 
images. With ALMA’s spatial resolution and cadence we may at last be able to resolve and 
identify the nature of these events and possibly answer some of the many open questions 
in this respect: Which physical processes generate the observed up- and down-flows and 
rotation in Explosive Events and which processes do produce these events in the first place? 
And do Explosive Events contribute directly or indirectly to the mass and energy transfer in 
the solar atmosphere? 


4.4 Solar flares 

Solar flare radio emission represents transient enhancements on top of a more slowly varying 
background (e.g. |Kundu|1965[|Kaufmann et al.|2004[ l. The great sensitivity and precision 
of ALMA means that we can study even the weakest events, and probe them throughout the 
corona and lower atmosphere. We discuss the possibilities in separate sections here for major 
events, the microflare-nanoflare domain, and the flaring lower atmosphere. See |Krucker et al.| 
( |2013| l for a recent review of observations and theory. 

4.4.1 Major events 

The microwave emission during flares comes primarily from mildly relativistic electrons via 
the gyrosynchrotron emission mechanism, with some role for thermal free-free emission. 
For a typical flaring loop filled with accelerated electrons, the standard gyrosynchrotron 
spectrum (see Sect. |2.4[ l reaches a peak around 5-20 GHz ( |Guidice & Castelli||1975[ |Nita| 
|et al.|[2004) l and decreases towards THz frequencies with a roughly power-law spectrum, 
whose spectral index is specified by the fast electron energy spectrum and their pitch-angle 
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anisotropy ( |Melrose||l968[|Ramaty||1969{ [Fleishman & Melnikov|[2003HKuznetsov et aT] 
[201 Ij l. In contrast, dense thermal plasma, which, for example, can be created as a result of 
chromospheric evaporation in response on the flare energy deposition, produces a flat flux- 
density continuum in its optically thin region and a Rayleigh-Jeans continuum at longer 
wavelengths. In some cases, available sub-THz observations do show either a decreasing 
spectrum, consistent with the non-thermal gyrosynchrotron mechanism, or an increasing 
continuum, consistent with the optically thick radiation. From this perspective, the discov- 


ery (jKaufmann et al.||2001[ Trottet et al.||2002t jLiithi et ^|2004[ [Kaufmann et al.||2004| 


jCristiani et al.|2008] Kaufmann et al.|2009a i of a continuum increasing with frequency and 

an apparently non-thermal component in the sub-THz range, turned out to be entirely unex¬ 
pected (as in Fig. left panel) even though entirely reasonable for a thermal one (jOhki &J 
|Hudson|1975| l. On a non-thermal hypothesis, it is currently unclear what particles (electrons, 
positrons, or maybe protons) are responsible for this spectral component, nor what emission 
mechanism is at work. Accordingly many ideas have been proposed recently ([ Fleishman &| 
|Kontar|2010[|Krucker et al.|2013[[Zaitsev et al.|2014l[Klopf et al.|2014l l. Nevertheless, no 
consensus has been reached concerning the relevant emission process yet since the currently 
available observations are made at only a few (in most cases, only two) widely separated 
frequencies with very limited spatial information. In order to distinguish the thermal and 
non-thermal emissions, the brightness temperature has to be measured across the whole 
source. The source sizes are estimated to be about 10-20” (see, e.g., [Silva et al.[199^ |R aulin| 
jet al.||T99^ [Ackermann et al.||2014[ [Kaufmann et al.[[2015[ for observations in different 
wavelength ranges) and thus are on the order of the ALMA primary beam width. Based on 
the existing data, including the limited infrared observations now available (e.g., [Xu et al.[ 
[2012[[Kaufmann et al. [2015) [Trottet et al.[2015[ l we expect to observe structure on all scales 
available to ALMA, and at all wavelengths. Interferometric observations with ALMA are 
therefore expected to greatly clarify the situation by adding the much-needed spectral and 
spatial information in the sub-THz range and thus enabling distinguishing between the ther¬ 
mal and non-thermal emission processes and the context association of the sub-THz source 
with other existing structures. As has been said, at least three distinct emission mechanisms 
are expected to make a contribution in the sub-THz emission from solar flares: thermal 
free-free, non-thermal gyrosynchrotron, and a non-thermal mysterious rising component. In 
principle, all these components may appear together, but presumably from distinct spatial 
locations, which calls for ALMA imaging of potentially complex sources. This complexity 
requires also the development of comprehensive modeling in three dimensions ( [Nita et JT] 
[20T5l l. 


Interestingly, the yet-to-be-explained distinct sub-THz component does not happen in 
each and every flare. It might not be present even in the most powerful events, like that 
shown in Fig. [T^ right. The emission is very strong and the spectrum peak frequency is 
very high (around 100 GHz), but no extra component is present in the sub-THz range. In 
such cases, ALMA will provide invaluable spectral information, which is the key to recov¬ 
ering the fast particle spectrum and other physical parameters of the flaring loops. [G^y et al.[ 
( [2013[ l employed a microwave spectrum in the range 1-18 GHz to develop a forward fit diag¬ 
nostics of the flaring loop parameters including the accelerated electron spectral slope. But 
they could neither recover the electron energy spectrum shape at the highest available energy 
nor determine the high-energy cut-off in that spectrum. However, this information is vitally 
needed to investigate the efficiency of the particle acceleration in flares and, thus, determine 
the responsible acceleration mechanism(s). At those high frequencies the gyrosynchrotron 
emission is produced by relativistic particles, which can be either accelerated electrons or 
secondary positrons (created in nuclear interactions). In order to distinguish the positron 
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Fig. 17 Illustration of the positron contribution effect on the polarisation of the gyrosynchrotron emission. 
Schematic energy spectra of fast electrons (red) and positrons (green) favourable for the high-energy positron 
detection are shown. The black line shows the composite spectrum dominated by electrons at low energies, 
but by positrons at high energies. The behaviour of the circular polarisation of the emission produced by 
these two components and by their composition is shown in the inset. A distinctive signature of the positron 
presence is the polarisation reversal at a high frequency shown by the red aiTow, which happens because 
the positrons preferably produce the ordinary mode emission in contrast to the electrons, which produce 
extraordinary mode emission. Figure from|Fleishman et al.H2013). 


and electron contributions, one has to measure the sense of circular polarisation along with 
direction of the magnetic field vector to determine the magneto-ionic mode dominating the 
emission-ordinary or extraordinary-to conclude in favor of positrons or electrons, respec¬ 
tively (see Fig.[T^ from |Fleishman et al.|2013[ l. 


So far, only indirect identifications of the wave mode are available with the Nobeyama 
data ( [Fleishman et al.|2013| l. ALMA with its broad spectral range, high spatial resolution, 
polarisation measurement purity (combined with the superior SDO vector magnetic mea¬ 
surements) will be an ideal tool for detecting the relativistic positron (and, thus, nuclear) 
component in solar flares — powerful, but almost unexplored diagnostics. In addition, after 
the positrons have become thermalized, they may form positronium, which might be de¬ 
tectable by ALMA. Thus, the evolution of the positrons from the original relativistic state 
to the final thermal state, when they interact with thermal electrons to form the positronium 
and eventually annihilate, can be nailed down. Then, detection of the thermal component 
from the flare can significantly constrain the energy balance, evolution, and partitions in 
flares. Although we cannot yet foresee what exactly will add the routine detection of the 
‘mysterious’ components, there is no doubt that ALMA observations of the flare emission 
will directly contribute to our understanding of the key questions about the solar flares. 
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4.4.2 Microflares and nanoflares 


The magnitudes of solar flares follow a power-law distribution over many decades (e.g. 
|Parnell & De Moortel|2012 and references therein). In this range the flare properties do not 
vary substantially (e.g. Qiu et al.|2004 |l, although there is a weak tendency for more energetic 
events to produce a harder (flatter) hard X-ray spectrum ( [Battaglia et al.|2005| l. In moderately 
active regions, microflares with non-thermal energies in excess of 10^'’ erg are found about 
once per hour (e.g., [Benz & Grigis|2002[[Krucker et al.|2002^ . At microwaves, non-thermal 
gyrosynchrotron emission has been suggested for a number of events (e.g. jFuerst et al.|1982[ 
[Gary et al.|1997[|Nindos et al.|199^|Kundu et al.|2006| l. 

Thermal flare signatures in the EUV down to a few times 10^^ ergs are found in quiet 
regions ( [Parnell & Jupp|[2000| l. [Krucker & Benzj l (2000l l reported frequent enhancements 
related to such quiet region flares in radio waves. Half of these events have a thermal radio 
spectrum (increasing with frequency). [Benz & Krucke'r| ( [1999[ l found thermal signatures of 
quiet region flares at a wavelength of 2 cm, originating mostly in the chromosphere. 

The great sensitivity of ALMA, even in single-dish operation, guarantees abundant ob¬ 
servations of these weaker events. At the bottom of the energy scale, with energy releases of 
order a few x 10“® those of major flares, we find the “nanoflare” concept. This interesting 
idea plausibly explains the entirety of coronal heating via tiny episodic flare-like processes 
([Parker|1988 i. At high occurrence rates the nanoflares would be individually undetectable 


Hudson 


1991[l, but indirect tests for their presence might reveal them; for lower rates 


(e.g. _ _ 

of occurrence one could possibly observe particular events (e.g. [Testa et al.|2014b) l. 


4.4.3 The lower atmosphere 


The presence of broad Her line profiles at the outer edges of expanding flare ribbons pro¬ 
vided early clues to the development of our current flare models, for example, and modem 
high-resolution observations reveal the presence of intense flare emission in white light in 
compact photospheric flare kernels ( [Krucker et al. [20 iTj [Martinez Oliveros et al.|2012 . An 
Her image of a C4.3-class flare, which occurred on 11 Nov 2010, is shown in Fig. 18 next 
to simultaneous images in the 17.1 nm and 30.4 nm passbands recorded with AIA/SDO. At 
ALMA frequencies one can observe the development of flare effects in the lower atmosphere 
with an almost entirely new capability. 

From the modeling point of view, flares or microflares offer a new challenge, but in prin¬ 
ciple a welcome one because of their “impulse response” excitation of physical processes. 
For example, [White et'af (1992 1 reported observations of an extraordinarily compact event 
(see also, e.g., Bastian et al.[2007[|Masuda et al.[2013| l, and [Jess et al. [ ( [20081 ) reported white- 
light emission from a GOES class Cl.2 flare on a 300km spatial scale. The IRIS mission 
has begun to provide imaging spectroscopy of the chromosphere, including of course flares 
(e.g. [Young et al.|2014] l, and thus has encouraged the development of extensive modeling 
tools as a routine part of the analysis of data. The scientific potential of simultaneous multi¬ 
wavelength observations of flares was already demonstrated during the Solar Maximum 
Mission (SMM) in the early 1980’s, during which several flares were observed simultane¬ 
ously in the UV at selected individual spectral lines, e.g. O V 137.1 nm. Si IV 139.4 nm etc, 
and at hard X-ray energies of several keV. Several papers reported simultaneous (to within a 
fraction of a second) increases in the UV lines and at hard X-ray energies (e.g., [Cheng et al.[ 
|1988[ l. It has now been shown ( [Doyle et al.[[2012] l that this is due to transient ionization, 
with the O V 137.1 nm line being enhanced in the first fraction of a second with the peak 
in the line contribution function occurring initially at a higher electron temperature than in 
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Fig. 18 Simultaneous imaging of a C4.3-class flare from 11 Nov 2010 at dilferent wavelengths. (Up¬ 
per left) Image taken in the wing of the Ha line (A A = 0.08 nm) at the Dunn Solar Telescope (NSO) with the 
IBIS instrument. (Upper right) AIA/SDO passband at 30.4 nm. (Bottom left) AIA/SDO passband at 17.1 nm. 
(Bottom right) A time-slice image in the Ho- showing the ribbon associated with the surge-like event, the ris¬ 
ing filament and its eruption as a dark feature. All images comprise 400 X 400 pixels with a pixel size of O.’T. 
The FOV of each image is therefore 40” X 40”. The images are based on the work by|Huang et al.|j2014a^. 


ionization equilibrium. The rise time and enhancement factor depend mostly on the elec¬ 
tron density. The high cadence capability of ALMA will be essential due to these transient 
ionisation effects. 


The appearance of strong impulsive-phase infrared sources ( [Xu et al.|2004 |at 1,560 nm 
and [Kaufmann et al.|2013| at 10,000 nm) establishes the major significance of the lower solar 
atmosphere in flare physics. These findings, together with the absence of evidence for strong 


beaming ( 

Kane et al.|1998HKontar & Brown|2006[|Martmez Oliveros et al.|2012HDickson| 

[& Kontar 

2013|l strongly suggest that ALMA will significantly contribute to answering im- 


portant open questions about the nature of solar flares. 


4.4.4 Quasi-periodicpulsations 


An intriguing and intensively debated feature of flaring energy releases are quasi-periodic 
pulsations (QPP) which are often seen in light curves of solar flares as well-pronounced 
periodic modulation of the emission intensity i Nakariakov & Melnikov|[2009| l. These pe¬ 
riodic variations were first detected by [Parks & Winckler ( 1969| l in hard X-rays and have 
been consequently observed in all wavelength bands including microwaves, white light, 
soft and hard X-rays, and gamma-rays. In a number of flares quasi-periodic pulsations are 
seen simultaneously with different instruments, ground-based and space-bome. For exam- 
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pie, quasi-periodic pulsations have been simultaneously detected in the microwave band 


1969 

Kane et al. 1983 [Nakajima et al. 

[19831 [Asai et al. [[20011 [Grechnev et al. 

2003b 

Fleishman et al.|2008 

Nakariakov et al.|2010 1, and in soft X-rays with broadband filters on 


Typical periods of quasi-periodic pulsations range from a fraction of a second to sev¬ 
eral minutes, and sometimes are of a very high quality. The modulation depth can reach 
100%. Statistical studies in the microwave band show that quasi-periodic pulsations are a 
common and perhaps intrinsic feature of flares ( [Kupriyanova et al.|2010[ l. The term “quasi”- 
periodic reflects that the variations do not always have a clear harmonic pattern, and often 
are amplitude- and period-modulated. This modulation can be connected with noise, or it 
can be associated with the gradual variation of the physical parameters of the flaring active 
region and hence carry useful information. 

It is commonly accepted that quasi-periodic pulsations can be produced by several non¬ 
exclusive mechanisms, such as modulation of the plasma parameters by MHD waves and 
oscillations (e.g. [Nakariakov & Melnikov| |2006|l, modulation of the non-thermal electron 
and positron kinematics by MHD oscillations (e.g. [Zaitsev & Stepan^ | 1982| l, periodic 
triggering of energy releases by external fast and slow magnetoacoustic waves ([Nakari-[ 
[akov et ar] [2006| [Chen & Priest[ 2006 1 , alternate currents (torsional Alfven oscillations) 
in current-carrying loops ( [Zaitsev et al. 1998[l, and self-oscillatory regimes of spontaneous 
magnetic reconnection (e.g. JKIiem^tal J2000] l, including the coalescence over-stability of 
two current-carrying loops ([Taj ima et al.[1987 1 . Quasi-periodic pulsations in the sub-second 
period range could be associated with wave-particle interactions ( [Aschwanden|1987 1 . Oc¬ 
casional observations of multi-period quasi-periodic pulsations (e.g. Inglis & Nakariak^ 
[2009[[^n Doorsselaere et al.|2011[ l might indicate that several different mechanisms could 
operate simultaneously, or that quasi-periodic pulsations are induced by several different 
co-existing MHD waves. Indeed, in some cases, multi-period quasi-periodic pulsations have 
been identified with standing slow and fast sausage modes of flaring coronal loops ( [Van[ 
[Doorsselaere et al.|2011[ |. 

The possible association of quasi-periodic pulsations with coronal MHD oscillations is 
based not only on the coincidence of the period ranges. In a number of cases quasi-periodic 
pulsations have the time signatures typical for the MHD waves detected in the corona in the 
EUV band. For example an exponentially decaying harmonic dependence was observed si¬ 
multaneously in the microwave quasi-periodic pulsations and in the EUV emission intensity 
oscillation caused by a standing slow magnetoacoustic wave ( [Kim et al.|2012[ l. However, 
in the majority of cases, it has been difficult to establish the specific mechanism respon¬ 
sible for a quasi-periodic pulsation because of insufficient spatial resolution. For example, 
the 10” resolution of the Nobeyama Radioheliograph (NoRH) does not allow resolving the 
width of a typical flaring loop and hence the spatial transverse structure of the quasi-periodic 
pulsations’s microwave source. ALMA’s high spatial resolution will certainly dramatically 
improve our understanding of this phenomenon. 

Observations of solar flares in the sub-THz band have shown high quality second and 
sub-second quasi-periodic pulsations (e.g [Makhmutov et al.|2003][Kaufmann et al.|2009a^ . 
Thus, revealing the nature of the dependence of the flaring emission intensity on the fre¬ 
quency in the GHz and sub-THz band would be incomplete without quantifying them in 
detail and explaining the nature of quasi-periodic pulsations in this band. 

Moreover, observations of quasi-periodic pulsations and their structure at different fre¬ 
quencies would provide us with additional important constraints necessary for the discrimi¬ 
nation between possible physical mechanisms responsible for the peculiarity of the sub-THz 
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spectram of flaring emission. One necessary step forward in this direction is the development 
and generalisation of forward modelling of microwave emission quasi-periodic pulsations 
performed by [Reznikova et al.| ( [2014| |2015[ l, including different sources of the periodic¬ 
ity, e.g. different modes of MHD oscillations, quasi-periodically varying parameters, such 
as the spectral, spatial and pitch angle distribution of the injected non-thermal particles. A 
preliminary analytical study of sub-THz quasi-periodic pulsations in terms of the alternate 
current model has recently been performed by |Zaitsev et al.| ( |2014| l and showed consistency 
with the available observations performed without the spectral and spatial resolution. In 
particular, this model could reproduce the observed dependence of the pulsation repetition 
rate on the signal intensity, which is different at the two separated frequencies of the obser¬ 
vation, 212 GHz and 405 GHz. The confident identification of the physical mechanism re¬ 
sponsible for quasi-periodic pulsations and, in the case of modulation by MHD oscillations, 
of the modulating MHD mode, requires information about the spatial and spectral struc¬ 
ture of quasi-periodic pulsations at different frequencies. For example, a decisive signature 
of the sausage oscillation are the anti-phase variations in different parts of the microwave 
band ( |Mossessian & Fleishman|2012| l. Also, the fine spectral structure of microwave quasi- 
periodic pulsations may contain unique signatures of different MHD modes ( |Yu et al.|2013^ . 
Forward modelling of similar time-spatial-spectral signatures of other possible mechanisms 
for quasi-periodic pulsations needs to be performed. 

Thus, high-precision imaging spectroscopy and polarimetry with ALMA with a spatial 
resolution of at least an order of magnitude better than NoRH’s and its high spectral resolu¬ 
tion, will revolutionise our understanding of flaring quasi-periodic pulsations. Furthermore, 
the quasi-periodic injection and acceleration of fast electrons (e.g. [Parks & Winckler|196^ 
|Kane et al.|1983t [Nakajima et al.|1983[ [Fleishman et al.|[2008[ l results in a corresponding 
measurable quasi-periodic impact of the electrons precipitating into the chromosphere (see 
Sect. [4.4.3| l. ALMA can measure the corresponding thermal chromospheric response and, 
thus, provide a much needed chromospheric diagnostics for the processes causing quasi- 
periodic pulsations. 


4.4.5 Particle beam heating of the chromosphere 


From the phenomenon of type III radio emission it is well known that particle beams form 
in the solar corona. Reverse drift bursts (e.g. |Aschwanden et al.|1995l l reveal electron beams 
moving downwards to the chromosphere, but it is still unclear how far down they go. Down¬ 
ward electron beams have a long theoretical development as part of the so-called thick-target 
model as an explanation for hard X-ray and white-light footpoint sources, and hence for flare 
energetics in general l [Kane & Donnelly 1 1971 [ [Brown| 1971 [ [Hudson| \912) . These electron 
beams are thought to hit the plasma at the loop footpoints in the lower (dense) solar atmo¬ 
sphere, where the bulk of hard X-ray (HXR) emission is produced. Recent observations with 
RHESSI have allowed to look into the structure of HXR footpoints (see, e.g., the review by 
Kontar et al. 20TTJ. The observations ( jAschwanden et al.|2002[[Kontar et al.[2010[[Battaglia[ 
& Kontar|2 011 1 established that higher energy HXR emission is indeed produced lower in 
the atmosphere. RHESSI observations ([Kontar et al.|2008[ [Saint-Hilaire et al.|20I0[ [Kon-[ 


[tar et al. [|20 1 Oj [Battaglia & Kontar|[20TT 2012 1 of hard X-ray foot-point sources suggest 

heights in the range 700 - 1200 km. Martinez Oliveros et al.| ( |2012[ l found significantly 
lower heights, at the highest (305 ± 170) km, but a further study of three (partially occulted) 
events ( [Kmcker et akpOlSl l suggests that the white light and HXR centroids are observed 
to be co-spatial and appear about 800 km above the photosphere. ALMA observations will 
provide crucial constraints on the formation heights. 
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Beam transport is also very attractive for explaining weaker events such as micro- and 
nanoflares (e.g., [Testa et al.|['2014b| l. It is commonly accepted that during solar flares par¬ 
ticle beams, which are accelerated by magnetic reconnection processes in the low corona, 
bombard the denser chromospheric layer. Owing to this bombardment the chromosphere is 
rapidly heated and the chromospheric plasma is evaporated into coronal parts of flare loops. 
Processes, which involve heating of the chromosphere by particle beams, have been stud¬ 
ied intensively by means of ID radiation hydrodynamic modeling ( jNagai & Emslie|1984 [ 
[MacNeice et al.|19^4t . 

Furthermore, it was also found that these processes are strongly influenced by the so- 
called return current ( [van den Oord|19^[Karlicky & Henoux|1992[ l, which is formed due 
to electromagnetic effects during the beam propagation through the flare atmosphere. The 
return current fully compensates the electric current of the beam and influences the location 
in the chromosphere, where the beam energy is deposited. It was proposed that the return 
current produces the impact polarisation of optical chromospheric lines ( [Karlicky & Henou^ 
[2002[ l. As detailed in Sect.j^ ALMA is expected to serve as a thermometer for the chromo¬ 
spheric plasma, which facilitates the study of all the processes connected with particle beam 
heating. However, ALMA alone will not be able to directly recognize particle beams but 
indirectly through changes in the emission due to beam heating of the chromosphere. In any 
case, ALMA observations need to be supported by further observations in a broad range of 
the electromagnetic spectrum including X-ray, gamma-ray and the THz range. 


4.4.6 Triggering mechanism of subflares in active regions 

In a recent work, [Tiwari et al.[ (pOT4| reported external triggering of solar subflares in a 
braided coronal magnetic structure observed by Hi-C ([Cirtain et al.|2013[ [Kobayashi et al.[ 
[2014[ l. Magnetic reconnection between two chromospheric/transition region loops below the 
braided coronal magnetic structure resulted in both a shorter and a longer loop. The shorter 
one presumably submerged, as implied by the flux cancellation seen at that site, while the 
longer loop erupted and interacted with the overlying braided coronal loop, thus triggering a 
subflare with an estimated thermal energy release of 10^® erg. At least ten such events were 
detected in four hours of observation using different channels of the Atmospheric Imaging 
Assembly (AIA) aboard the Solar Dynamics Observatory (SDO). This observation suggests 
a new possible mechanism for the external triggering of subflares in addition to the one 
most commonly thought, namely the spontaneous internal triggering as a consequence of 
excessive energy accumulation by footpoint shuffling of magnetic field lines ( [Parker|1988[ (. 
It remains to be seen, however, how common this external triggering is and how impor¬ 
tant it is for coronal heating in general. High resolution chromospheric images obtained by 
ALMA, in particular when co-observing with future Hi-C flights, would provide important 
constraints on the magnetic and thermal topology of the active region chromosphere before, 
during and after subflares, and the external triggering of these events. 


4.5 Chromospheric heating in regions with strong magnetic field 

Maintaining the high temperatures in the magnetic chromosphere, i.e. above network and 
plage regions, requires a heating rate that is larger than that of the corona ( [Withbroe &[ 
[Noyes|1977[ ), yet the dominant mechanism for heating the magnetic chromosphere remains 
unknown. A variety of observations have revealed evidence for potential candidate heating 
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mechanisms fro m spicules (|de Pontieu e t al.|2007{|Rutte n|2006[l, Alfven waves (|De Pon-| 


tieu et al.|2007[pess et al.| 2009[ Wedemeyer-Bohm et al.|2012{ De Pontieu et a l.|2014a^ 

’ ^ ([Jefferii ' ' 


magneto-acous tic shocks ([Jefferies^ et al.||2006[|Hansteen et al.|2006||Beck et al.||2008[ l to 

gravity waves ( [Straus et al.||2008 1. Similarly, a variety of theoretical models for chromo¬ 
spheric heating in magnetic regions have been studied, including high frequency waves 
([Hasan & van Ballegooijen|2008|l, Alfven waves (van Ballegooijen & Cranmer|2010[|Song| 


([Hasan 

[& Vasyli unas|201 4|l, transve rse kink waves ( Van Doorsselaere et al.|2008| l, multi-fluid ef¬ 
fects ( Krasnoselskikh et al.||2010j [Martmez-Sykora et al.||2012| l, and plasma instabilities 
(e.g., |Fontenla et al.|2008 1. Waves and oscillations and related implications for atmospheric 
heating are discussed in more detail in Sect. |4.6| 

This large number of possible sources of heating has made it difficult to determine which 
of these many processes actually dominate(s). This is in part because many of the mod¬ 
els lack the necessary physical complexity, including multi-fluid effects, non-LTE radiative 
transfer, or time dependent ionization. However, equally importantly the requirements on the 
heating (i.e., the “heating rate”) have often been derived from simplified, semi-empirical ID 
models (e.g., |Fontenla et al.|1993| l that differ significantly from the time-dependent, highly 
spatially structured solar chromosphere. In other words, a model’s “success” has often been 
determined by comparing the energy dissipation from a mechanism to an average heating 
rate derived from a ID static atmosphere. 

To make further progress it is critical that more quantitative and realistic comparisons are 
made between observations and models, e.g., through comparisons of observations (e.g., in 
the Mg II and Ca II lines) and synthetic observables from models (|de la Cruz Rodriguez et al.| 
|2013a[[Leenaarts et al.|2013j|Pereira et al.|2013| l. While many of the chromospheric spectral 
lines require advanced calculations that include non-LTE radiative transfer, the availability 
of ALMA observations (which largely avoid non-LTE calculations) will provide a critical 
new component to constrain the spatio-temporal properties of chromospheric heating in and 
around magnetic network and plage regions for a wide range of heights. 

ALMA observations have the potential to provide more direct constraints for theoret¬ 
ical modelers that, when combined with sub-arcsecond resolution observations of chro¬ 
mospheric lines observed from the ground or space (e.g., with IRIS and/or SST/CRISP/ 
CHROMIS), will lead to new insights concerning which of the many forms of mechanical 
energy that have been invoked by models actually dominate the deposition of non-radiative 
energy in the magnetized chromosphere. Because coronal loops are rooted in the network 
and plage regions, a better understanding of chromospheric heating has direct implications 
for coronal heating models. 


4.6 Chromospheric oscillations and waves 
4.6.1 Wave propagation in the solar atmosphere 

Waves and oscillations are interesting not only from the point of view that they can propagate 
energy into the chromosphere and dissipate that energy to produce non-radiative heating, but 
they also carry information about the structure of the atmosphere in which they propagate. 
Since the sound speed is a slowly varying function of temperature in the chromosphere, 
pressure disturbances with frequencies above the acoustic cut-off propagating upwards in 
the chromosphere provide information about the difference of the formation height of the 
accompanying temperature and velocity perturbations, or, reversely, if the difference of for- 
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mation heights is known, allow inferences about the phase speed of the waves and hence 
about the wave mode. 

Several authors jM ein|1971]|Mein & Mein|1976[|Lites et al.|1982t[Staiger et al.|1984[ 


[Fleck & Deubner||1989 1 who studied the wave propagation behavior in the chromosphere 

by analyzing time series of the Call infrared triplet lines at 854.2nm and 849.8 nm found 
zero phase differences between the Doppler shifts of the two lines. The vanishing phase 
lag between the two lines was interpreted as evidence of a standing wave pattern in the 
chromosphere caused by total reflection of the waves at the transition region ( |Mein|197T| 
|Mein & Mein|[T976l [Staiger et al.||198^ [Fleck & Deubner]|1989| l. This conclusion was 
based on the assumption that the Doppler signal of the two Ca II infrared lines is formed 
at different heights in the chromosphere, an assumption that was later demonstrated to be 
wrong ( [Skartlien et al.|1994[ l. [Fleck et al.[([1994a[l, [Fleck et al.[p994b^ and [Hofmann et ah] 
( [1995[ l extended the earlier work of Fleck & Deubner ( 1989 1 and included the CaIIK and 
He 11083 nm lines together with the Ca II 854.2 nm line in their observations. It is generally 
believed that the He 1 1083 nm line is formed high up in the chromosphere just below the 


transition region (e.g., Brajsa et al.||l996 2007 Jurdana-Sepic et ^|2009 


and references 

therein) and that the formation height of Ca IIK 3 also is different from that of the Ca II in¬ 
frared triplet ( [Skartlien et al.|1994^ . The phase difference spectra between these three lines, 
however, still showed a nearly uniform value close to 0 °, especially at frequencies above the 
cutoff frequency, providing further evidence of a dominating non-propagating component of 
the chromospheric wave field. 

Quasi-simultaneous observations of diagnostics with different formation heights, as will 
be possible in a linear manner with ALMA, have been carried out by [Doyle et al.[ ( [r999[ l, 
who observed the two lower-chromospheric lines NI 131.9nm and CII 133.5nm with the 
SUMER instrument on the SoHO spacecraft. The differences between (intensity and veloc¬ 
ity) power spectra of network and internetwork regions are in line with results based on 
Call power spectra, whereas it seems that CII intensity lags that of NI by up to 16 sec near 
3 mHz when phase spectra are averaged along the slit. 

According to the dispersion relation of acoustic-gravity waves in a stratified atmosphere, 
propagating sound waves display a linear increase of the phase difference with frequency. 
Interestingly, the analysis of recent IRIS observations of the Mg II h and k line also shows 
vanishing phase differences between the ks, k 2 v and k 2 r CallK line core features, which are 
believed to form at different heights. Why is the apparent phase speed of high frequency 
acoustic waves so high? Are these results misleading because of complex radiation transfer 
effects in these optically thick lines? 

[Loukitcheva et ah] ( [2004[ [2006[ [2008[ l, and [White et al.[ ( f2006| l have demonstrated the 
feasibility of measuring chromospheric oscillations in the mm range. Multi-wavelength time 
series of ALMA observations of the temperature fluctuations of inter-network oscillations 
should allow travel time measurements between different heights as these disturbances prop¬ 
agate through the chromosphere and thus should finally settle the long-standing question 
about the propagation characteristics of acoustic waves in the chromosphere. 


4.6.2 Alfven waves 


In addition to acoustic waves, Alfven waves may play an important role in chromospheric 
heating. In the chromosphere, collisions between ions and neutrals (see also Sect. |4.2.5| are 
an efficient dissipative mechanism for Alfven waves with frequencies near the ion-neutral 


collision frequency (De Pontieu et al.[[2001 Khodachenko et al.|2004 

Leake et aL|2005[ 

[Zaqarashvili et al.[|2011 

Russell & Fletcherl 2013| Soler et al. [2013b 

ai. Estimations of 
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the heating rate due to Alfven waves damped by ion-neutral collisions ( [Song & Vasyliun^ 
[201 l[|Goodman|201 l[|Tu & Song|2013| l suggest that this mechanism may generate sufficient 
heat to compensate the radiative losses at low altitudes in the solar atmosphere. Recently, 
jSoler et ah] (| 2015| ( investigated the damping of Alfven waves as a function of height in a 
static chromospheric model based on the semi-empirical model F of jFontenla et al.| ( |l993| (. 
They found a critical interval of wavelengths for which impulsively excited Alfven waves 
are overdamped as a result of the strong ion-neutral dissipation. Equivalently, for period¬ 
ically driven Alfven waves there is an optimal frequency for which the damping is most 
effective. The critical wavelengths and optimal frequencies vary with height and are very 
sensitive to the ion-neutral collision cross section for momentum transfer used in the theo¬ 
retical computations. Uncertainties in the collision frequency and cross section cause uncer¬ 
tainties in the various transport coefficients that govern basic collisional phenomena in the 
plasma ( jMartmez-Sykora et al.|201^ . On the other hand, jZaqarashvili et al.| ( |2013| ( inves¬ 
tigated the role of neutral helium and stratification on the propagation of torsional Alfven 
waves through the chromosphere. While high-frequency waves are efficiently damped by 
ion-neutral collisions, jZaqarashvili et al!] ( |2013[ l concluded that low-frequency waves may 
not reach the transition region because they may become evanescent at lower heights ow¬ 
ing to gravitational stratification. Hence, propagation of Alfven waves through the chromo¬ 
sphere into the solar corona should be considered with caution. The ability of ALMA to 
observe different layers in the chromosphere by using different bands, along with the high 
temporal and spatial resolutions, may be crucial to understand how Alfven waves actually 
propagate and damp in the chromosphere and thus how they contribute to chromospheric 
heating. ALMA offers the unprecedented opportunity to study high-frequency waves in the 
chromosphere. 

Magnetohydrodynamic waves are also ubiquitous in solar prominences ([Oliver & Ballesterj 
[2002[[Lm et al.[2009t[Hillier et al.[2013^ . Since the properties of the prominence plasma are 
akin to those in the chromosphere, dissipation of wave energy due to ion-neutral collisions 
may also play a role in the energy balance of prominences. In the same way as in the chro¬ 
mosphere, ALMA observations may shed light on the impact of high-frequency waves on 
prominence plasmas. 


4.6.3 Resonant absorption and associated heating 


Transverse MHD waves have been shown to be ubiquitous in the solar atmosphere ([DeJ 
[Pontieu et al. [2007] [McIntosh et al.[201 l[[Okamoto et al.|2007[[Tomczyk et al.|2007) l. Such 

waves constitute important heating agents due to their relatively high energy flux. Numerous 
theoretical and numerical studies have shown that resonant absorption is an expected phe¬ 
nomenon accompanying the propagation of transverse MHD waves along chromospheric or 
coronal waveguides ( lonsoiJ2978|^ Sakurai et al.|199f [Arregui et al.|2008[ [Terradas et al.[ 


|2008b[[Van Doorsselaere et al.|2008 Verth et aL|2010[ i. An observational consequence of 

this mechanism is the damping of transverse motions of the waveguides (such as flux tubes) 
in the plane-of-the-sky (POS) due to the conversion of energy between the kink mode and 
Alfven waves. This conversion occurs at the resonance points, where the phase speed of 
the kink wave is equal to that of the Alfven waves, and may not involve significant dissipa¬ 
tion. ft has been therefore difficult to associate the observed wave damping with heating and 
to estimate the exact energy contribution from these waves to chromospheric and coronal 
heating. 

High-resolution ALMA observations at chromospheric and coronal heights (including 
prominences and coronal rain) can contribute in several ways towards narrowing down the 
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heating contribution of transverse MHD waves. For instance, the high spatial resolution 
of ALMA allows to refine the upper limits for the energy contents of these waves, available 
prior to dissipation or mode conversion ( |Goossens et al.|20T3j|Van Doorsselaere et al.|2014^ . 
Another important way is by clarifying the dissipation mechanisms of these Alfvenic waves, 
which must occur at presently unresolved spatial scales. As explained earlier, resonant ab¬ 
sorption can lead to strong damping of the POS transverse motion and a corresponding 
increase of the azimuthal flows in the boundary layer, which behave like m = \ torsional 
Alfven modes. Through coordinated observations between ALMA and spectrometers such 
as IRIS, both motions associated with transverse MHD waves can be detected (POS and 
azimuthal motions along the line-of-sight), allowing precise MHD mode identification and 
correct application of MHD seismology techniques (|Verth et al.||2011| ?;|Okamoto et al.| 
|2015b| l. 

Also, resonant absorption is expected to generate small spatial scales in which, through 
the process of phase mixing, heating is produced. Furthermore, analytical and numerical 
modelling demonstrates that low amplitude transverse MHD modes can have a significant 
impact in the morphology of coronal or prominence flux tubes l|Ofman et al.|1994{ |Te rradas| 
|et al.|2008a}|Soler et al.|20r0l l. The strong velocity shear obtained through resonant absorp¬ 
tion can lead to Kelvin-Helmholtz instabilities (KHI), which extract the energy from the 
resonance layer and dump it into heat through viscous and ohmic processes in the generated 
vortices and current sheets. The Kelvin-Helmholtz instability vortices that result show-up 
as strand-like or thread-like structures in coronal or prominence flux tubes, respectively 
l |Antolin et al.|[20T4l |2015a[ [Okamoto et al.|[?015b| l. These vortices should be detectable 
at a resolution of 20 % of the flux tube’s radius, which would be ~ 0.3" for a flux tube 
with 1000 km radius and thus an achievable goal with ALMA. Detection of such Kelvin- 
Helmholtz instability vortices and the ensuing turbulence would allow to correctly estimate 
the wave energy being dissipated. 


4.7 Magnetic loops in the upper atmosphere 
4.7.1 The fine-structure of coronal loops 


For a long time solar coronal loops have been considered to consist of bundles of thin strands 
that remain undetected with the current instrumental capabilities (see for e.g., [Gomez et al.| 
|1993[|Reale|2010] and references therein). Today, that perspective continues to remain as 
prevalent as ever. Coronal loops were first detected in coronagraphic observations in the 
1940s ( [Bray et al.|1991[|Aschwanden|2004l and references therein). Within coronal loops, 
neighbouring field lines are considered to be thermally isolated and substructures are re¬ 
ferred to as “strands”. It is likely that most observations represent superpositions of hun¬ 
dreds of unresolved strands at various stages of heating and cooling ( |Klimchuk|2006| l. Other 
studies suggest that elementary loop components should be even finer, with typical cross- 
sections of the strands on the order of 10-100 km ( [Beveridge et al.|2003[ |Cargill & Klimchukj 
|2004[ [DeForestj 2007 1 . See also the results obtained from imaging spectroscopy at radio 
wavelengths by |Chen et al.| ( [2013[ l with the Karl G. Jansky Very Large Array (VLA), which 
support loop diameters of less than 100 km. In order to fully describe the sub-structures of 
coronal loops one fundamental issue needs to be addressed, namely the unknown spatial 
scale of the coronal heating mechanism l [Reale|2010[ (. The multi-stranded interpretation of 
a coronal loop cross-section predicts a broad temperature distribution of unresolved threads 
([Aschwanden & Boerner|2011^. [Peter et al.|([2013[l concluded that the fundamental reason 
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656.3 (+0.1032) nm 



Fig. 19 Left: A image of a post-flare loop taken with the CRisp Spectro- 

Polarimeter (CRISP) at the Swedish 1-m Solar Telescope (SST) in the Ho- far red-wing (656.3 nm 
+0.1032nm) (Scullion et al.|20T4) . The appearance of coronal rain (dark strands in Her) is co-spatial with 
high temperature emission in the SDO/AIA 17.1 nm channel (yellow contours, log T = 5.8) from the coronal 
loop. Right. Histogram of coronal rain widths for two datasets combining various wavelengths: Her with the 
SST (in red, average measured logarithmic temperature log T x 3.8), CallH with Hinode/SOT (orange, 
log T X 4), Mg Ilk 279.6nm with IRIS/SJI (light green, log T » 4.2), CII 133.0nm with IRIS/SJI (dark 
green, log T k 4.3), Si IV 140.0 nm with IRIS/SII (blue, log T k 4.8), He II 30.4 nm with SDO/AIA (light 
purple, log T « 5), FeIX 171 with SDO/AIA (dark purple, log T » 5.9), and FeXII 19.3 nm with SDO/AIA 
(black, log T « 6.2). Explanation for these datasets can be found in |Antolin et al.|f2015b( . The bin size in 
each wavelength is set to the spatial resolution of each instrument. Strand-like structure in these loops extends 
from chromospheric to at least transition region temperatures. The shape of the distributions strongly suggest 
that the bulk of the coronal rain widths remains undetected at present day resolution. Also compare to the 
statistical analysis of strands in this coronal loop by|Scullion et al.|j2014). 


for the constant cross-section is due to the presence of a small temperature gradient across 
the structure, hence, a simpler spatial structure would suffice. [Schmelz et al.|([201 ij l studied 
cool loops that were chosen in the 17.1 nm channel from SDO/AIA and found 12 loops with 
narrow temperature distributions, consistent with isothermal plasma that could be modelled 
independently as a monolithic loop. In 2012 the High-resolution Coronal Imager (Hi-C: 
ICirtain et al.|2013l l imaged the 1.5 MK solar corona at the highest ever resolution (in that 
passband) of 216-288 km (Winebarger et al.]|2014^ . From this data, [Brooks et al!] ( |2013| l 
measured the Gaussian widths of 91 coronal loops and the resulting distribution had a peak 
at a cross-sectional width of 270 km. 

Coronal loops in active regions can easily be subject to the thermal instability mecha¬ 
nism, through which the plasma undergoes catastrophic cooling from coronal temperatures 
down to chromospheric temperatures ( jAntiochos et al.|1999t|Aschwanden et al.|2001| l. The 
end result of this process is coronal rain, chromospheric material observable by ALMA. As 
we will see in the next sections, coronal rain provides unique advantages for solar atmo¬ 
spheric research, such as being a high resolution tracer of the local and global magnetic 
field topology of loops. 

4.7.2 Coronal rain - a proxy for the spatial and thermal inhomogeneity of the solar corona 

A major question in coronal heating concerns the spatial distribution of the heating mech¬ 
anism in loops. There is substantial evidence for most of the heating in active region loops 
being concentrated at the footpoints, which implies that thermal conduction is transporting 
the energy from energetic events at the footpoint upwards along the loops and is thus the 
main mechanism maintaining the corona in active regions l|Klimchu^|2006[ |Reale|[20T0). 
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Numerical studies have shown that such spatial distribution of heating often sets the loops 
into a state of thermal non-equilibrium (Antiochos & Klimchuk 1^91^ Antiochos^etaUl^iQ^ 


IMiiller et al.|2003l|2004|[2005b[|Mendoza-Briceno et al.|2005 |Antolin et al.|2010HSusino| 


|et al.|20im|Mikic et al.|2013l l, in which cycles of heating (where loops become over-dense) 
are followed by cooling (in which loops deplete) due to the thermal instability mechanism 
( |Parker| 195 3 [ |Field| 1965 [ I . The lowest temperatures in these cycles depend on a broad pa¬ 
rameter space, and range from cool chromospheric temperatures to transition region temper¬ 
atures. During the cooling stage the plasma rapidly condenses into cold and dense clumps 
termed coronal rain, often observed falling along loop-like trajectories in a large number of 


Brekke[[1998 

Schrijver[[2001| 

Muller et al. |2005a Ugarte-Urra et al.[[2006[ O’Shea et al. 

2007[ Tripathi 

et al.|2009| Antolin & Rouppe van der Voort|2012| Antolin et al.[2015b I. 


Recent studies suggest that coronal rain is a fairly common phenomenon of active re¬ 
gions playing an important role in the chromosphere-corona mass cycle ( |Antolin & Roupp^ 
[van der Voort|201^|Liu et al.|2012[|Antolin et al.|2015b[ l. While these studies suggest that 
the temperatures in coronal rain can go down to 2000 K or so, it is not yet clear whether 
such low temperatures are common in the rain. Through the thermostat and high resolu¬ 
tion capability of ALMA, the atmosphere above an active region can be probed for thermal 
bremsstrahlung from the cool material. Also, coronal rain is expected to have strong emis¬ 
sion in recombination lines that are detectable with ALMA. Determining the presence of 
such cold material in the solar atmosphere is of large importance since it would imply that 
coronal heating is an extremely inhomogeneous process, with a very cold and frequent coun¬ 
terpart. 

Combining both a chromospheric and coronal nature, coronal rain constitutes the highest 
resolution window into the substructure of loops. At the highest spatial resolution achiev¬ 
able today, the coronal rain clumps appear organised in strand-like structures with ever de¬ 
creasing widths on the order of a few 100 km or less (see Fig. [T^ right panel), suggest¬ 
ing a tip-of-the-iceberg scenario in which loop substructure may not be fully resolved yet 
(Antolin et al.|20f^|Antolin & Rouppe van der Voort|2012[ [Scullion et al.|2014]|Antolin] 


et al.|2015b i. These results have been supported by numerical simulations of coronal ram 

formation, through which a higher percentage of coronal rain clumps are obtained at sub¬ 
resolution scales ( [Fang et ar]|2013[ l. The strand-like structure of coronal rain is not only 
limited to chromospheric temperatures but extends to at least transition region temperatures 
( [Antolin et al.|2015bl l. Such results point to the existence of fundamental substructure within 
coronal loops, consistent with the theoretically predicted strands (see Sect. |4.7.T| i. By pro¬ 
viding high spatial resolution maps of sizes and temperature of coronal rain, ALMA can 
help determining the bulk of the size distribution of such fundamental substructure and its 
role in the chromosphere-corona mass cycle. 

The nature of coronal rain strands and prominence threads remains elusive ([Tandberg-J 
|Hanssen|1995[|Lin|201 l[|Engvold|2015| see also Sect. |4.8| l. The formation and evolution of 
such structure may be strongly linked to coronal heating mechanisms. Recent work suggests 
the generation of such fine-scale structure and associated heating through dynamic instabil¬ 
ities (such as KHI [Antolin et al.||2014| [2015a[ [Okamoto et al.p015b[ l. With high spatial 
and temporal resolution, and multi-passband observations, ALMA can clarify the genera¬ 
tion of such fine-scale structure by tracking the heating events in these cool structures and 
by determining the thermal and non-thermal processes involved. 

Perhaps the most famous example of coronal rain is post-flare loops, in which the strong 
basal heating from the flare produces large chromospheric evaporation leading to over-dense 
strongly radiating loops which end up cooling catastrophically and falling back to the chro- 
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mosphere ( |Foukal|1976[|1978] l. This multi-thermal aspect of the plasma is often invoked to 
explain the usually large flux densities observed in the millimetre and sub-millimetre range 
l |Trottet et al.|201 i'l l but is not well constrained during the flare event. Also, the spatial dis¬ 
tribution of the main thermal component during a flare is not well known. ALMA can help 
determining the amount and spatial distribution of the chromospheric temperature plasma 
such as coronal rain contributing to the thermal radiation. Observations of coronal rain with 
ALMA will help determining the thermal and spatial inhomogeneity of the solar corona, and 
thus establish fundamental aspects of coronal heating and thermal instability mechanisms. 


4.8 Prominences and filaments 


As already mentioned in Sect. 4.1 prominences (aka filamentj^ are cool, dense regions 
of plasma lying in a much hotter and rarified coronal environment, which are supported by 
coronal magnetic fields against gravity. Comprehensive reviews of prominence physics can 
be found in |Tandberg-Hanssen| i 1995[ l, Labrosse et al. (|2010|l, |Mackay et al.| (|2010|l, |Parenti| 
l |2014[ l or in the recent book by | Vial & Engvoldp 2015 I. 

The primary contribution of ALMA to the understanding of solar prominences is its 
unique capability to measure the temperature of the cool prominence plasma with high spa¬ 
tial and temporal resolution (see, e.g., |Loukitcheva et al^ |2004t [Heinzel & Avrett| |2012[ 
|Perez-Le6n et al.||2013 I. Previous studies based on mm and sub-mm radio emission (see, 


e.g., Vrsnak et al.|1992 [Bastian et al.|1993^|Harrison et al.|19^|Irimajiri et al.|1995^ had 

to rely on low-resolution observations available at the time. Now ALMA will be able to re¬ 
solve the thermal fine-structure of prominences on far smaller spatial scales, perhaps reach¬ 
ing the physical dimensions of the prominence fine structures. Moreover, the high temporal 
resolution of ALMA will facilitate studying the evolution and dynamics of the (thermal) 
the small-scale structures, including, e.g., localized plasma flows, oscillations and passing 
waves. 

Basic plasma parameters can easily be obtained, assuming that the prominence plasma 
is optically thin, that the background signal can be safely removed, and that the radiometric 
calibration is correct ( |Perez-Le6n et al.|2013[ l. As the thermal structure of the cool promi¬ 
nence plasma is still not well understood, it will be of primary interest to probe the plasma 
under different viewing angles (e.g. filament on the disk vs. prominence above the limb), 
and across the wavelength range offered by ALMA. Due to the increase of opacity of the 
prominence plasma and thus the variation of the formation region of the continuum with 
increasing wavelength (see, e.g., [Heinzel et al.||20f5l l, ALMA will be able to effectively 
scan the prominence thermal structure from the coolest fine-structure cores towards the hot¬ 
ter layers surrounding them. This will provide new insights into the structure of the lower 
temperature prominence-corona transition region (PCTR), which surrounds a whole promi¬ 
nence or individual prominence fine structures, and resulting in direct implications for the 
energy balance problem in prominences. 

ALMA observations of the thermal structure of prominences can be tested and cross- 
calibrated using synthetic temperature maps obtained by forward modelling. For this pur¬ 
pose, multi-dimensional prominence models with realistic temperature distribution includ¬ 
ing the PCTR can be used, e.g., the 2D models of |Gouttebroze & Labrosse| (|2009 1 and 
[Heinzel & Anzer| (|2001|l, or the 3D whole-prominence fine structure models of Gunar &| 


Prominences are refen'ed to as (dark) filaments when observed in absorption against the (bright) solar 

disk. 
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Fig. 20 Left: Modelled (ideal) brightness temperature map at 100 GHz. The white circle shows ALMA field- 
of-view for single pointing, the crosses are the centres of pointing for interferometric mosaic. Fourty-nine 
pointings are necessary for Nyquist sampling of the displayed part of the prominence. Right: Corresponding 
image as seen by ALMA at 3 mm in Band 3 simulated with the CASA::simobserve() task. Images from 
|Heinzel et al.HTOfS) . We note that the dynamic range of images obtained by interferometric observing at the 
limb may be challenging. 


|Mackay| ( |2015[ l. [Heinzel et al.| ( |2015| l simulated the visibility of solar prominences and the 
detectability of their internal structures in future ALMA observations. Unlike in the case 
of the quiet solar chromosphere (see Sect. |4.2[ l, the brightness temperature model is based 
on the relation between the kinetic plasma temperature and the derived integral intensity of 
the Ha line (see |Heinzel et al.| 2015| and references therein). Figure|^shows the modeled 
(i.e. ideal) brightness temperature map calculated from relations (10) and (14) in |Heinz^ 
|et al.| ( |20T5l l together with pointings for an ALMA interferometric mosaic at 3 mm (ALMA 
Band 3) and the corresponding image reconstructed from the ALMA visibilities. 


According to previous radio observations of filaments l |Raoult et al.|1979] l, the intensity 
contrast with respect to the quiet sun decreases with frequency, being ~5 % at 3.5 mm (see 
also |Chiuderi Drago et al.|2001^ . At higher frequencies, [Bastian et al.| ( [1993a[ l found negli¬ 
gible contrast at and near the location of Ha filaments at 0.85 mm. This finding demonstrates 
that filaments can expected to be observable in the low frequency ALMA bands, with the 
TP system being necessary due to their large extent. As a matter of fact, a filament corridor 
is visible in the SW quadrant in the ALMA map shown in Fig.[^ In addition, dense patches 
of filamentary material might be observable at higher frequencies. Such observations will 
provide information on the density and temperature structure and, possibly, on the magnetic 
field of filaments, complementary to the information obtainable from prominences. 

In summary, the unique capabilities of ALMA will result in significant contributions to 
our understanding of solar prominences. For instance, the anticipated results concern (i) de¬ 
tailed information on the thermodynamical parameters the prominence fine structure plasma, 
(ii) the nature of the PCTR and its role in conducting energy from the corona to the promi¬ 
nence, (iii) the differences in 3D magnetic structures between prominences, and (iv) how 
the prominence plasma collects and evolves (see also the associated issue of coronal rain in 
Sect. |4.7.^. Fur thermore, ALMA will help advance the understanding of wave propagation 
(cf. Sect. 4.6. l| l and damping in prominences and with it more details of the prominence 
plasma and magnetic field properties. ALMA will also help to investigate currently debated 
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Fig. 21 Characteristic sizes and timescales of different rotating magnetic phenomena in all scales: giant 
tornadoes (squares), rotating network fields (RNFs: diamonds), magnetic tornadoes (triangles, tip downward), 
explosive events (EEs: circles), spicules (triangles, tip upward), and an extremely large tornado (ELT) on solar 
atmosphere and, in addition, astrophysical jets in young stellar objects (YSOs: bowties), in X-ray binaries 
(XRBs: sandglasses), and in active galactic nuclei (AGNs: stars). The individual data points for objects in the 
solar atmosphere are explained in jWedemeyer et al.|2013^ and for objects in YSOs, XRBs, and AGNs are 
explained in the text in Section [4.9| Filled boxes indicate the areas of observed data points whereas unfilled 
boxes indicate the area of theoretical possibility. 


‘giant’ solar tornadoes (e.g., Orozco Suarez et al.|2012[[Su et al.|2012[|Li et al.|2012[|Pane-| 


|sar et al.|2013[[Wedemeyer et al.|2013a[ and references therein), which are often associated 

with solar prominences. While recent observations allowed more information to be gathered 
on the state of the plasma l [Levens et al.|2015j , open questions remain as to the exact nature 
of the magnetic structure supporting this kind of phenomenon. 


4.9 Potential implications for other fields of astrophysics 

The Sun serves as a natural plasma laboratory, which hosts a multitude of physical processes 
that can also occur in other astrophysical objects either under similar conditions or on other 
scales. Many of the results that are expected from ALMA observations of the Sun, therefore 
have potential implications for other helds of astrophysics. 

Observations of the (sub-)millimeter emission from the solar atmosphere can be com¬ 
pared to (spatially unresolved) ALMA observations of other stars, thus opening up a much 
larger parameter space. Given its high sensitivity, ALMA is well-suited for such stellar ob¬ 
servations as demonstrated by, e.g.,|Liseau et al.|((2015|l. They managed to resolve the stellar 
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components of the binary a Centauri with ALMA and measured the flux density in three 
different bands, thus providing important constraints for the thermal structure of the atmo¬ 
spheres of these two sun-like stars. Comparing such stellar ALMA observations to solar re¬ 
sults promises important clues for the heating and activity of stellar atmospheres in general, 
including, e.g., the source of the so-called “basal flux” (e.g., |Schrijver|19?7}|Perez Martmez| 
|et al.|20lT| and flares (see Sect. |4.4[ l. For example, quasi-periodic pulsations (see Sec. 4.4.4^ 


are found both for solar and stellar flares. With dedicated campaigns ( |Anfinogentov et ^ 
2013 I and with the Kepler mission ( [Walkowicz et al.|2011[[Balona|2012[ Shibayama et al. 


2013 I, flares and quasi-periodic pulsations in them are regularly detected. However, seis¬ 


mology with these pulsations is even less straightforward, since the spatial scales of the 
flare region cannot be directly measured on the stars. The correct identification of the mech¬ 
anism responsible for quasi-periodic pulsations in solar flares will enable a more targeted 
seismology of the stellar flares as well. As in the solar case, ALMA will also be an excel¬ 
lent instrument for the study of stellar flare oscillations. For instance, [Welsh et al.|p006l l 
detected oscillations of 30 to 40 sec in flares observed with the GALEX satellite in several 
stars. These oscillations were interpreted as due to magneto-hydrodynamic waves in active 
region coronal loops (cf., e.g., [Van Doorsselaere et al.|2011[ |. 

Observations of magnetic tornadoes and vortex flows in the upper solar atmosphere (see 
Sect. |4.2.4| i seem to suggest a relation between the characteristic size and the lifetime of 
different kinds of rotating magnetic field structures on the Sun (Fig.[2T|l. This finding gives 
rise to the question if such a relation could be extrapolated to larger scales on the Sun and 
elsewhere in the universe and if there is a physical mechanism which would explain such 
a scalability. Large-scale, long-lived rotating magnetic fields are indeed observed as part of 


Bacciotti et al.|2002| Coffey et al.|2004 

Carrasco-Gonzalez et al.|2010|l, through galactic 

compact objects (Mirabel & Rodriguez 1998 Fomalont et al. 2001 [Gallo et al.[2005|[Begel- 

man et al. 2006b to active galactic nuclei i 

Junor et al.|1999 Lobanov & Zensus|2001[ Beck- 

ert & Falcke|2002| 

Gebhardt & Thomas 

2009 Mahmud et al.|2013 i and a galactic centre 

molecular tornado i 

Sofue|2007b (see Fig. 21 1 . Although the basic physical mechanism that 


creates and drives these jets is yet to be understood, rotating magnetic fields play a critical 
role in creating, driving, and directing jets across all scales. Actually, successful theoretical 
models, which have been developed for these different systems, incorporate magnetic fields 
threading the central objects such as accretion discs or rapidly spinning stars and compact 
objects. These models are tested observationally by measuring polarised radio emission (see 
[Pudritz et al.|2012l for a review and references therein) and compared to laboratory exper¬ 
iments (for an overview see |Bellan et al.|2009| and references therein). ALMA will enable 
magnetic field measurements as well as vortex velocities, which are indispensable for un¬ 
derstanding this type of phenomena and for testing if the suggested scalability and potential 
similarities between the individual phenomena outlined above are true or not. 


5 Conclusions and Outlook 

While ALMA has already produced some remarkable results for a wide range of astrophys- 
ical topics, its potential for solar observations is yet to be revealed. Observing the Sun re¬ 
quires the development of observing modes, which are quite different from the usual modes. 
In return, such observations will produce some fascinating results. Essential contributions 
are anticipated for the science cases outlined in this review article. In particular, impor¬ 
tant steps towards advancing our understanding for central questions in solar physics are 
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expected, for instance, concerning the heating of the outer solar atmosphere, solar flares, 
prominences and their implications for stellar activity and space weather. 

It is important to remember that ALMA observations will be even more valuable when 
put into context with other observations of the chromosphere, e.g. obtained with IRIS (Mg II 
lines), SDO (magnetograms, continuum images, He II line), other ground-based solar obser¬ 
vations (Htr, Ca IIH and K lines. He 11083 nm, magnetic field maps, G-band, etc), as well as 
forthcoming observatories (e.g. Solar-C). Taken all together, new advances can be expected 
to address fundamental questions about the Sun. We expect solar observing proposals to be 
accepted beginning with ALMA cycle 4 in 2016. ALMA is expected to remain operational 
for several decades, and, as most other ground-based observatories, will further improve 
its already impressive capabilities with time, thus extending the range of possible scientific 
applications in the future. 
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